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ABSTRACT 


A  five  tube  column  test  condenser  was  modified  to  facili¬ 
tate  easier  tube  removal  and  installation  and  to  allow  for 
the  simulation  of  larger  depth  tube  bundles.  Heat  transfer 
coefficients  were  then  determined  for  sixty-two  runs  conducted 
at  steam  supply  pressures  of  35,  40,  45,  and  5G  psig  with  the 
system  operated  at  test  condenser  pressures  of  approximately 
2  and  15  psia.  Each  run  was  five  minutes  long  with  tempera¬ 
ture,  pressure,  and  flow  rate  data  taken  at  one  minute  intervals. 

The  results  revealed  that  the  average  tube  wall  tempera¬ 
ture  and  the  heat  transfer  coefficients  for  a  given  tube  of  a 
tube  bundle  increased  as  the  mean  vapor  velocity  increased, 
but  decreased  as  the  ar. ount  of  condensate  inundation  increased. 
The  experimentally  obtained  values  for  the  heat  transfer  coef¬ 
ficients  were  compatible  with  the  Nusselt  predictions  for  the 
runs  conducted  at  15  psia,  but  greatly  exceeded  the  Nusselt 
predictions  for  the  runs  conducted  at  2  psia.  However,  for 
the  2  psia  runs,  the  heat  transfer  coefficients  for  the  top 
tubes  agreed  closely  with  an  experimental  correlation  proposed 
by  Fujii.  The  simulation  of  a  tube  bundle  of  more  than  five 
active  tubes  per  col'imn  at  the  2  psia  test  condenser  pressure 
demonstrated  heat  transfer  coefficients  closer  to  the  Nusselt 
theory  than  when  only  five  active  tubes  were  considered. 

Recommendations  to  improve  the  test  apparatus  and  to  con¬ 
duct  additional  tests  are  piovided. 
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NOMENCLATURE 


C 

p 

d 

d 

o 

g 

h 

h. 

1 


hM 

Mu 


h 


o 


2 

Heat  transfer  area  of  one  tube  (m  ). 

2 

Heat  transfer  area  of  i  th  tube  (m  ). 

o 

Mean  flow  area  of  condenser  cross-section  (m“). 
Specific  heat  at  constant  pressure  (kj/kg  •  °C) . 
Specific  heat  at  constant  pressure,  i  th  tube 
(kJ/kg  •  °C) . 

Specific  heat  at  constant  pressure  (kJ/kg  •  C) . 
Outside  diameter  (m) . 

Outside  diameter  (m) . 

Outside  diameter  (m). 

2 

Acceleration  of  gravity  (9.81  m/s  ). 

?  o 

Heat  transfer  coefficient  (kw/m“  •  C) . 

Experimentally  determined  value  for  the  heat  transfer 
coefficient  of  the  i  th  tube  (kw/m2  •  °C). 

Latent  heat  of  vaporization  (kj/kg) . 

Local  heat  transfer  coefficient  for  the  N  fn 
tube  ( kW/m2  •  °C) . 

Average  heat  transfer  coefficient  for  a  column 
of  N  tubes  (kW/m2  •  °C)  . 

Heat  transfer  coefficient  calculated  from  the 
Nusselt  equation  (kW/m1-  •  C). 

Heat  transfer  coefficient  derived  from  the  Fujii 
correlation  (kW/m2  •  °C)  . 
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k 


k 


L 


m 


rai 


m 


CQND 


M 


M 


Mu 


m 


Fr 


D 

‘  TC 
q 

q/A 


ReL 

S/D 


Thermal  conductivity  (w/m  •  °C) . 

Thermal  conductivity,  liquid  (w/m  •  °C). 

Mass  flow  rate  of  cooling  water  (kg/s). 

Mass  flow  rate  of  cooling  water  for  the  i  th  tube 
(kg/s)  . 

Total  mass  flow  rate  of  the  condensate  (kg/s). 

The  number  of  tubes  in  a  column,  or  the  tube  number 
of  a  given  tube . 

The  number  of  unit  ceils  across  the  width  of  a 
condenser . 

Mean  Musselt  number  for  a  tube. 

Mean  Musselt  number  for  a  tube  without  inundation. 
Prandtl  number. 

Supply  steam  pressure  (psig) . 

Condenser  pitch  to  diameter  ratio. 

Test  condenser  pressure  (mm  Hg  abs.). 

Heat  transferred  per  unit  time  (kW) . 

Heat  transferred  per  unit  time  for  the  i  th  tube  (kW) . 
Heat  flux  (kW/m^). 

Two  phase  Reynolds  number,  dQ/v^. 

Spacing  to  diameter  ratio,  same  as  Prp. 

Cooling  water  inlet  temperature  (°C) . 

Cooling  water  outlet  temperature  (°C) . 

Saturation  temperature  of  steam  (°C) . 

Vapor  (steam)  temperature  (°C)  . 
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w 


Average  tube  wall  temperature  (°C). 

UM  Vapor  velocity  (m/s). 

3 

v  Specific  volume  of  the  steam  (m  /kg). 

Volumetric  flow  rate  of  water  to  the  desuperheater 
( ml /min) . 

V  Vapor  velocity  based  on  the  mean  flow  area  (m/s). 

Vp^,  Volumetric  flow  rate  of  water  to  the  porous  tube  (m/s)  . 

V$C  Volumetric  flow  rate  of  condensate  from  the  secondary 

condenser  (ml/min) . 

V_  Volumetric  flow  rate  of  condensate  from  the  test 

i  W 

condenser,  condensation  rate  (ml/min). 
w  The  amount  of  condensate  formed  on  the  N  th  tube  per 

unit  time  as  used  in  Equation  (5);  the  rate  of  inunda¬ 
tion  falling  onto  a  tube  as  used  in  Equation  (10). 

W  The  amount  o.f  condensate  draining  onto  the  N  th  tube. 

Equation  ( 5 )  . 

w°  The  rate  of  condensation  on  a  tube  corresponding  to 

Nu°,  Equation  (10). 

AT  Temperature  difference  (°C)  . 

C  Acceleration  parameter  (kAT/yh^-g) . 

2 

v  Kinematic  viscosity  (m  /hr). 

U  Dynamic  viscosity  (kg/m-s). 

S  Heat  capacity  parameter  (CpT/h^) 

p  Density  (kg/m^) . 

3 

Pv  Vapor  density  (kg/m  ) . 
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I.  HISTORICAL  BACKGROUND 

A  continuing  interest  exists  in  the  Navy  to  design  smaller 
and  lighter  marine  steam  condensers  by  the  use  of  enhanced 
heat  transfer  techniques.  This  is  an  outgrowth  of  recent  im¬ 
provements  in  turbine  machinery  and  boiler  design  which  have 
brought  about  increases  in  the  horsepower  to  weight  ratio  of 
marine  propulsion  plants,  but  with  the  lack  of  similar  improve¬ 
ments  made  to  condenser  design.  Through  the  application  of 
enhanced  heat  transfer  techniques  to  condenser  design,  further 
significant  increases  in  the  horsepower  to  weight  ratio  can  be 
attained  and/or  steam  cycle  efficiency  can  be  further  increased. 

Marine  steam  condenser  design  is  currerrtly  based  upon  the 
Heat  Exchange  Institute  (HEI)  standards  for  steam  condensers 
[Ref.  l]  and  the  standards  of  the  Tubular  Exchanger  Manufacturers 
Association  (TEMA)  [Ref.  23.  Search  [Ref.  3]  studied  existing 
marine  steam  condensers  to  determine  the  feasibility  of  improv¬ 
ing  condenser  performance  by  enhanced  heat  transfer  techniques 
and  concluded  that  these  techniques  could  result  in  a  forty 
percent  reduction  in  condenser  weight  and  volume,  indicating 
a  need  to  establish  new  design  criteria  for  marine  steam 
condensers . 

Beck  [Ref.  4],  Pence  [Ref.  5],  Reilly  [Ref.  6],  Fenner 
[Ref.  7],  and  Ciftci  [Ref.  3]  conducted  experimental  research 
into  various  kinds  of  enhancement  using  a  single  tube  heat 
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transfer  apparatus .  The  above  investigations  determined 
that  for  the  same  diameter  tube,  the  overall  heat  transfer 
coefficients  of  enhanced  tubes  exceeded  those  for  smooth 
tubes  by  almost  100%. 

Marto,  Reilly,  and  Fenner  [Ref.  9]  determined  that  most 
of  the  improvement  occurred  on  the  cooling  water  side  of  the 
tube  due  to  a  combination  of  increased  surface  area  and  an 
increased  turbulence  and  swirl  in  the  cooling  water  flow. 

This  gain  in  heat  transfer  was  accomplished  at  the  expense 
of  a  substantial  increase  in  the  cooling  water  pressure  drop 
due  to  friction  factors  as  much  as  ten  times  larger  than  for 
smooth  tubes .  They  also  noted  that  enhanced  tubes  with  deep 
grooves  had  the  best  heat  transfer  performance,  and  that  for 
a  given  tube  geometry  at  constant  groove  depth,  heat  transfer 
performance  was  dependent  upon  the  groove  pitch  or  helix  angle. 

To  fully  evaluate  the  outside  heat  transfer  performance 
of  enhanced  tubes,  it  was  required  to  consider  the  effects 
of  condensate  inundation  and  vapor  velocity.  Demirel  [Ref.  10] 
and  Eshleman  [Ref.  Ill  consecutively  modified  the  test  appara¬ 
tus  originally  designed  by  3eck  (for  single  tube  experimenta¬ 
tion)  to  study  the  effects  of  condensate  inundation  on  a 
column  of  five  horizontal,  smooth  tubes  in  a  staggered  arrange¬ 
ment.  Their  work  was  inconclusive  however,  since  problems 
with  non-condensable  gases,  side  drainage,  and  secondary  steam 
flow  within  the  condenser  resulted  'n  ambiguous  tube  perfor¬ 
mance  data  for  the  bottom  two  tubes  . 
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A  new  test  apparatus  was  designed  by  Morrison  [Ref.  12] 
in  an  attempt  to  eliminate  the  problems  encountered  by  Demirel 
and  Eshleman,  and  to  provide  a  system  which  would  operate 
closer  to  the  conditions  of  the  Nusselt  theory.  Morrison’s 
test  condenser  was  smaller  in  length  than  that  of  Demirel  and 
Eshleman  and  employed  a  new  diffuser  section  at  the  steam  inlet 
in  an  attempt  to  eliminate  the  secondary  steam  flow.  The  con¬ 
denser  tubes  were  also  arranged  in-line  vice  staggered  to 
assist  in  the  reduction  of  side  drainage.  The  above  factors 
also  allowed  for  more  efficient  removal  of  non-condensable 
gases  by  the  steam  flow  through  the  condenser.  Thermocouples 
were  also  installed  within  the  active  tube  walls  so  that  direct 
calculation  of  the  outside  heat  transfer  coefficients  was 
possible.  Preliminary  data  taken  after  the  publication  of 
Morrison's  work  indicated  that  the  apparatus  was  suitable  for 
condensate  inundation  studies,  producing  tube  performances  in 
close  accordance  with  the  Nusselt  theory. 

Though  Morrison's  test  apparatus  eliminated  or  reduced  the 
problems  encountered  by  Demirel  and  Eshleman,  several  other 
problems  were  noted.  Air  leaks  at  the  joints  and  welds  limited 
the  vacuum  that  could  be  attained  and  introduced  non-condensable 
gases  into  the  system.  The  tube  wall  thermocouples  were  very 
fragile  and  broke  easily  when  tubes  were  installed  or  removed. 
This  was  compounded  by  the  fact  that  the  0-ring  tube  sealing 
method  used  required  considerable  force  to  be  applied  to  the 
tubes  during  installation  and  removal.  Dropwise  condensation 
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persisted  on  the  upper  two  tubes  despite  the  use  oi  various 
tube  cleaning  methods,  and  the  difficulty  of  tube  removal 
and  reinstallation  limited  the  frequency  of  thorough  tube 
cleaning. 

The  objectives  of  this  thesis  were  therefore  to: 

1.  Further  modify  Morrison's  test  apparatus  to  eliminate 
or  reduce  the  tube  wall  thermocouple  breakage  and  the 
active  tube  removal  and  reinstallation  problems. 

2 .  Establish  a  tube  cleaning  procedure  to  eliminate  or 
reduce  the  dropwise  condensation  problem. 

3.  Further  modify  Morrison's  test  apparatus  through  the 
inclusion  of  a  porous,  water  supply  tube  above  the 
first  of  five  active  tubes  such  that  the  apparatus 
could  be  used  to  simulate  tubes  in  larger  tube  bundles. 

4.  Obtain  baseline  heat  transfer  performance  data  for 
the  test  condenser  utilizing  16  urn  O.D.  smooth  copper 
tubes . 
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II.  THEORETICAL  BACKGROUND 


Nusselt's  theory  of  film  condensation,  based  on  work  com¬ 
pleted  in  1916,  provided  the  basis  for  most  of  the  knowledge 
developed  to  date  in  regard  to  condensation  on  plates  and 
horizontal  tubes.  Nusselt's  work  was  based  on  the  following 
assumptions  for  single  tubes  as  noted  by  Nobbs  [Ref.  13]; 

1.  The  wall  temperature  is  constant. 

2.  The  flow  is  laminar  in  the  condensate  film. 

3.  Heat  transfer  in  the  condensate  is  by  conduction, 
and  subcooling  may  be  neglected. 

4 .  The  fluid  properties  are  constant  within  the  condensate 
film. 

5.  The  forces  due  to  hydrostatic  pressure,  surface  tension, 
inertia,  and  vapor-liquid  interfacial  shear  are  negli¬ 
gible  when  compared  to  the  viscous  and  gravitational 
forces . 

6 .  The  surrounding  steam  and  vapor/liquid  interface  are 
at  saturation  temperature. 

7.  The  film  thickness  is  small  when  compared  with  normal 
tube  diameters  and  the  effects  of  curvature  are  small. 

When  a  column  of  tubes  is  to  be  considered,  Eissenberg 
[Ref.  14]  noted  that  the  following  two  added  assumptions  are 
required : 
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8 .  Condensate  drains  as  a  laminar  sheet  from  a  tube  bottom 
to  a  tube  top  such  that  velocity  and  temperature  gra¬ 
dients  are  not  lost  in  the  fall  between  tubes. 

9 .  The  saturation  temperature  and  the  tube  wall  tempera¬ 
ture  are  constant  for  all  tubes  in  the  bank. 

For  a  single  horizontal  tube,  Nusselt  determine'3’  that  th<= 
outside  heat  transfer  coefficient  could  be  obtained  from  the 
equation : 


hM 

Nu 


.72  5 


K3  P(  P  -  ».,)  hfg  g  1/4 

M  D  (T  -  T  )  -* 

s  w 


(1) 


where , 

k  =  thermal  conductivity  of  the  condensate 

p  =  density  of  the  condensate 

Pv  =  density  of  the  vapor 

hfg  =  latent  heat  of  vaporization 

g  =  acceleration  of  gravity 

U  =  dynamic  viscosity  of  the  condensate 

D  =  outside  diameter  of  the  tube 

Ts  =  temperature  of  the  steam 

=  average  tube  wall  temperature 

Jacob  [Ref.  15]  determined  that  for  a  tube  column  of  N 
tubes,  an  average  outside  heat  transfer  coefficient  for  the 
column,  h^ ,  could  be  obtained  from  the  equation: 
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Using  equations  (1)  and  (2),  the  Nusselt  result  is  normally 
expressed  as  the  ratio  of  the  average  heat  transfer  coefficient 
(h^)  to  the  heat  transfer  coefficient  of  the  top  tube  (h^  ): 


!h-  =  (r1/u 

hNu 


(  3) 


It  may  also  be  expressed  as  the  ratio  of  the  local  heat  transfer 
coefficient  of  a  given  tube  in  a  column  to  the  heat  transfer 
coefficient  of  the  top  tube: 


N3/4  -  (N  -  i)3/1+ 


(4) 


It  has  been  noted  that  most  experimentally  determined 
results  are  on  the  order  of  twenty  percent  greater  than  those 
obtained  from  the  above  equations.  This  is  due  to  the  nature 
of  the  assumptions  made  in  the  derivation  of  Nusselt ' s  theory. 
In  actual  marine  steam  condensers,  the  vapor  can  have  a  fairly 
high  velocity  over  a  large  portion  of  its  flow  path;  and  under 
changing  turbine  speeds,  the  steam  velocity  is  certainly  not 
negligible.  The  vapor  velocity  causes  vapor-liquid  interfa¬ 
cial  shear.  With  a  downward  flow  of  vapor  and  condensate,  a 
shear  force  is  added  to  the  gravitational  force .  Consequently 


20 


as  the  condensate  film  velocity  increases,  the  film  thickness 
decreases  and  the  heat  transfer  coefficient  from  the  vapor 
to  the  wall  increases. 

Moreover,  as  noted  by  Kutateladze  [Ref.  16];  Shklover  and 
3uevich  [Ref.  17];  Yung,  Lorenz,  and  Ganic  [Ref.  18];  and 
Eissenberg  [Ref.  14],  condensate  does  not  flow  from  each  tube 
as  a  continuous  film  or  sheet.  Instead,  the  condensate  flows 
in  the  form  of  separate  drops  or  streams.  The  condensate 
dripping  from  an  upper  tube  will  split  around  the  tube  below 
but  will  not  flow  axially;  it  may  hit  the  lower  tube  off- 
center  or  may  miss  the  lower  tube  altogether,  depending  upon 
the  tube  bundle  geometry.  The  thickness  of  the  film  on  a 
lower  tube  can  therefore  vary  greatly,  and  where  thinnest, 
greatest  heat  transfer  will  occur. 

Fuks  [Ref.  19]  developed  an  experimental  correlation  for 
the  normalized  local  heat  transfer  coefficient  in  the  form: 


h,.  . .  -0.07  n  m 

»  =  c  :  =  >r0-07  (5) 

hi  « 

Where  W  is  the  amount  of  condensate  draining  onto  the  N  th 
tube  per  unit  time  and  w  is  the  amount  of  condensate  formed 
on  the  N  th  tube  per  unit  time. 

Chen  [Ref.  20]  derived  a  formula  which  considered  boundary 
conditions  influenced  by: 


21 


1.  the  momentum  gain  of  the  falling  condensate  between 
tubes ,  and 

2.  the  condensation  of  vapor  on  the  condensate  between 
tubes . 

His  result  was: 


-  1/4 

it  :  M--/4  [l  ♦  0.2CCM-1)  [  Mttt  ]1  ( 6 ) 

nNu  *•  1*0. 95  5-0.15;; 


where  , 

5  =  KAT/uhf 
5  =  CpAT/h 

and , 


fg 


5  =  5/Pr. 

Eissenberg  [Ref.  14  ]  experimentally  investigated  the  effects 
of  steam  velocity,  condensate  inundation,  and  non-condensable 
gases  on  the  heat  transfer  coefficient.  Based  on  the  results 
of  those  experiments,  he  developed  a  side  drainage  model  for 
staggered  tube  bundles  by  which  condensate  may  drain  onto  tubes 
not  directly  below  one  another  as  depicted  in  Figure  1.  For 
comparison  to  the  above,  Figure  2a  depicts  the  Nusselt  in-line 
tube  bundle  drainage  model  in  which  condensate  drains  to  lower 
tubes  in  sheets;  while  Figure  2b  illustrates  a  more  realistic 
in-line  tube  bundle  drainage  model  in  which  condensate  drains 
to  lower  tubes  in  drops,  with  considerable  splashing. 

Eissenberg  identified  the  following  conditions  as  affecting 
side  drainage: 
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1.  Orientation  -  Staggered  tube  bundles  are  more  suscep¬ 
tible  to  side  drainage  than  in-line  tube  bundles; 

2.  Spacing  -  The  smaller  S/D  is,  the  more  frequently  side 
drainage  should  occur; 

3.  Momentum  -  The  greater  the  horizontal  component  of 
momentum  of  a  drop  leaving  a  tube,  the  greater  the 
side  drainage  ; 

4.  Steam  velocity  -  When  steam  flows  horizontally  across 
tubes  at  sufficient  velocity,  the  drop  trajectory  will 
reflect  the  added  lateral  momentum.  When  steam  flows 
vertically,  its  direction  change  with  each  tube  may 
also  impart  lateral  momentum  to  the  drops;  and 

5.  Misaligned  tubes  -  A  tube  misaligned  in  a  bundle  may 
receive  greater  or  lesser  amounts  of  condensate  depend¬ 
ing  on  its  orientation  with  respect  to  the  side  tubes. 

Utilizing  a  pure  side  drainage  model  for  a  staggered  tube 
bundle,  Eissenberg  obtained: 

M  1  /  U 

r-—  =  0.69  +  0.42  N  x/  (7) 

hMu 


Since  experimental  results  deviate  from  each  other  and 
from  the  Musselt  theory,  results  are  often  fitted  to  the 
following  equation  form: 


h 


N 


N 


-s 


(8) 


2  3 


■-*"**► 


where  s  is  an  experimentally  determined  exponent  with  a 
reported  range  of  values  of  0.07  to  3.20. 

Nobbs  [Ref.  13]  concluded  the  following  from  his  investi¬ 
gation  of  the  combined  effects  of  downward  vapor  velocity  and 
condensate  inundation  on  the  rate  of  condensate  formation  in 
horizontal  tube  bundles: 

1.  Vapor  velocity  increases  the  condensate  heat  transfer 
coefficient  on  both  inundated  and  uninundated  tubes 
in  a  tube  bundle  . 

2.  The  effect  of  inundation  is  to  reduce  the  heat  transfer 
coefficient.  The  rate  of  reduction  with  increase  in 
inundation  rate  becomes  smaller  as  the  vapor  velocity 
is  increased. 

3.  The  condensate  drainage  path  is  often  not  vertically 
downwards  but  in  a  diagonal  direction.  This  can  result 
in  tubes  receiving  differing  amounts  of  inundation. 

Fujii  [Ref.  21]  used  the  data  obtained  in  the  experiments 
of  Nobbs  and  Mayhew  to  correlate  the  effects  of  inundation  and 
vapor  velocity.  For  the  top  or  uninundated  tube  in  an  in-line 
tube  bundle,  the  following  equation  was  obtained: 

Mu°  =  10.74  Re.0,312  (9) 

m  L 

where 

Nu°  =  Nusselt  number,  h  d  /kT  ,  for  steam  without 
m  ’  o  o  L  ’ 

inundation 
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ReT  =  two-ohase  Reynolds  number,  U  d  /v, 

k^  =  thermal  conductivity  of  liquid 

d  =  outside  diameter  of  tube 
o 

=  kinematic  viscosity  of  liquid 
U  =  vapor  velocity 

For  any  tube  in  an  in-line  tube  bundle,  Fujii  also  obtained 
the  equation: 

0.65 


where  , 

Mu  =  mean  Nusselt  number  for  a  tube 
m 

w  =  rate  of  inundation  falling  onto  a  tube 

w°  =  rate  of  condensation  of  a  tube  corresponding  to 


When  calculating  vapor  velocity,  most  studies  have  used 
the  entrance  vapor  velocity  which  is  based  upon  the  condenser 
cross-sectional  area.  However,  as  steam  flows  through  a  con¬ 
denser,  its  velocity  is  alternately  increased  and  decreased 
due  to  the  restriction  in  actual  flow  area  caused  by  the  tubes. 
Therefore,  using  the  entrance  velocity  (or  for  that  matter, 
the  maximum  velocity  obtained  at  the  point  of  minimum  flow 
area)  does  not  result  in  a  true  picture  of  the  effect  of 
vapor  velocity  on  heat  transfer.  Nobbs  [Ref.  13]  used  a  mean 
velocity  of  vapor  based  upon  a  mean  flow  width.  Figures  3 
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and  4  show  the  derivations  of  the  mean  flow  width  for  the 

cases  of  a  square  pitch  and  an  equilateral  triangle  tube 

bundle,  respectively.  From  these  Figures  one  can  obtain  for 

the  in-line  or  square  pitch  tube  bundle: 

A  _  =  N  dL  (PT  -  t/4P_)  (11) 

mt  c  l  T 

wnere , 

A  c  =  mean  flow  area 
mf 

N  =  number  of  unit  cells 

c 

d  =  tube  outside  diameter 

P^.  =  tube  pitch  to  diameter  ratio 

L  =  tube  length 

Then,  for  the  purposes  of  mean  vapor  velocity  calculations 
and  for  use  in  Equation  (9),  the  following  was  used  in  this 
study : 

v  =  m-C0MD-.-y  (12) 

m  A  f 
mf 

where  , 

=  mean  vapor  velocity 
mC0ND  =  mass  rate  of  condensate 

v  =  specific  volume  of  vapor  at  the  steam  saturation 
tempe  rature 

Note  that  although  the  above  is  an  improvement  over  the  use 
of  the  entrance  vapor  velocity,  it  does  not  account  for  the 
decrease  in  vapor  velocity  through  the  tube  bundle  due  to  the 
reduction  in  the  amount  of  steam  flow  due  to  the  condensation 
of  the  steam. 
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Figure  5  compares  the  equations  of  Nusseit  (Equation  3), 
Eissenberg  (Equation  7),  and  Chen  (Equation  6)  for  fi^/h^. 
Figure  6  compares  the  equations  of  .'Jusselt  (Equation  4)  and 
Fuks  (Equation  5)  for  h^/h^ . 
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Ill . 


EXPERIMENTAL  APPARATUS 


The  test  apparatus  designed  by  Morrison  [Ref.  13]  was 
modified  to  allow  for  the  simulation  of  an  active  tube  column 
having  a  number  of  tubes  in  increments  of  five  tubes  deep 
(five,  ten,  fifteen,  etc.  tubes  deep)  and  to  correct  the  resign 
deficiencies  noted  above.  Some  components  of  the  original  test 
apparatus  were  also  relocated  to  allow  for  more  efficient  opera 
tion  of  the  system.  The  following  system  descriptions  update 
those  presented  in  Morrison's  thesis  to  incorporate  the  modifi¬ 
cations  that  were  made.  Figures  7a  and  7b,  and  the  key  thereto 
show  the  assembled  test  apparatus  as  modified  and  detail  the 
major  components. 

A.  TEST  CONDENSER 

The  test  condenser,  shown  in  Figures  8-12,  measured  305  mm 
x  305  mm  x  79  mm  and  was  made  of  stainless  steel.  The  dimen¬ 
sions  allowed  for  a  maximum  of  twenty-seven  16  mm  O.D.  tubes 
arranged  in  an  in-line  configuration  of  three  columns  of  nine 
tubes  each.  However,  the  in-line  tube  configuration  initially 
used  was  one  of  a  column  of  five  active  tubes  flanked  on 
either  side  by  a  column  of  five  dummy  tubes.  Directly  above 
the  upper  a-cive  tube,  a  porous,  distilled  water  supply  tube 
was  positioned,  flanked  on  either  side  by  dummy  tubes. 
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The  square  in-line  arrangement  of  tubes  with  a  spacing-to- 
diameter  ratio  of  1.5  was  utilized  for  the  conducted  experiments. 
However,  the  test  section  design  allowed  for  the  rearrar.gment 
of  dummy  tubes  to  give  a  staggered  arrangement  as  well.  This 
was  made  possible  by  machining  slots  in  the  end  plates  of  the 
test  condenser  and  providing  tube  sheets  for  the  dummy  tubes 
which  fit  into  the  slots  such  that  they  could  be  changed  to 
suit  the  type  of  arrangement  desired.  Each  end  of  the  dummy 
tubes  was  inserted  into  its  respective  tube  sheet,  and  the 
tube  sheets  were  then  slid  into  the  end  plate  slots  and  secured 
into  place  by  screws. 

An  opening  (Fig.  12)  running  vertically  along  the  center- 
line  of  each  condenser  end  plate  provided  for  active  and  porous 
tube  penetration  of  the  test  section.  The  tubes  were  then 
fixed  in  position,  with  a  spacing-to-diameter  ratio  of  1.5, 
by  means  of  nylon  tube  sheets  that  were  attached  to  the  exter¬ 
ior  of  the  condenser  end  plates.  A  new  set  of  tube  sheets 
would  be  required  fur  each  desired  spacing-to-diameter  ratio. 

To  facilitate  ease  of  active  tube  removal  and  installation, 
and  to  prevent  damage  to  the  thermocouples  soldered  into  the 
active  tube  walls,  the  tube  sheets  and  tube  sealing  method 
were  redesigned.  For  each  test  condenser  side,  a  nylon  tube 
sheet  one  inch  thick  was  fabricated  with  six  holes  (S/D=1.5) 
such  that  about  0.5  mm  tube  clearance  existed  which  allowed 
for  the  tubes  to  be  easily  slid  in  or  out  of  the  tube  sheets. 
Grooves  were  machined  into  the  exterior  side  of  each  tube  sheet 


around  each  hole  to  accommodate  O-rings .  The  aluminum  tube 
sheets  from  the  original  test  condenser,  but  with  a  sixth  hole 
added  and  all  holes  machined  to  about  a  0 . 5  mm  tube  clearance, 
were  then  used  as  sealing  plates.  When  installed,  rubber 
gaskets  were  placed  between  the  test  condenser  end  plates  and 
the  nylon  tube  sheets  ;  0-rings  were  placed  in  the  tube  sheet 
grooves;  and  the  aluminum  seal  plates  were  positioned  on  the 
tube  sheets.  These  assemblies  were  then  secured  to  the  test 
condenser  end  plates  by  means  of  screws.  With  the  screws 
loosened,  the  active  tubes  could  be  readily  removed  or  in¬ 
stalled;  but  when  the  screws  were  tightened  down  with  the 
tubes  in  place,  the  0-rings  were  sufficiently  compressed 
around  the  tubes  to  provide  a  tight  seal. 

The  diffuser,  exhaust  plenum,  transition  piece,  vortex 
annihilator,  and  the  exhaust  and  condensate  piping  (Figs. 

3-10)  as  well  as  the  viewing  window  were  unchanged  from 
Morrison's  design. 

3.  TEST  CONDENSER  TUBES 

The  active  tubes  were  made  of  16  mm  O.D.  half  hard,  smooth 
copper  tubes  with  a  1.65  mm  wall  thickness.  In  Morrison's 
original  apparatus,  the  active  tubes  extended  through  the  tube 
sheets  such  that  the  portion  of  the  tubes  extending  outside 
of  the  test  condenser  acted  as  highly  conductive  fins.  New 
active  tubes  were  made  approximately  305  mm  long,  the  width 
of  the  test  condenser,  with  16  mm  O.D.  stainless  steel  tube 
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extensions  fixed  to  each  end  of  the  active  tubes  for  the  pur¬ 
pose  of  tube  sealing  and  connection  to  the  cooling  water  supply 
and  discharge  lines. 

Four  thermocouples  were  soldered  into  grooves  90  degrees 
apart  within  the  wall  of  each  active  tube  such  that  the  sensing 
portion  of  each  thermocouple  was  about  at  the  midpoint  of  the 
wall  thickness  .  The  grooves  extended  from  the  center  of  each 
copper  section  axially  about  15  cm  into  the  stainless  steel 
extension.  Upon  completion  of  the  soldering,  the  soldered 
surfaces  of  the  tubes  were  ground  and  tubes  polished  so  that 
the  tube  surfaces  were  as  uniformly  smooth  as  possible.  The 
high  thermal  conductivity  of  the  copper  tubes  and  silver  solder 
ensured  that  the  thermocouples  gave  reasonably  accurate  measure 
ment  of  the  tube  wall  temperature.  The  free  ends  of  the  thermo 
couples  were  secured  to  the  tubes  with  tape  to  minimize  the 
risk  of  damage. 

The  dummy  tubes  flanking  the  active  tubes  were  made  of  16 
mm  O.D.  stainless  steel.  These  tubes  served  to  direct  the 
steam  flow  so  as  to  simulate  actual  conditions  in  a  condenser. 
Cooling  water  was  not  supplied  to  these  tubes  and  they  did  not 
penetrate  the  test  condenser  end  plates . 

C.  POROUS  TUBE  WATER  SUPPLY  SYSTEM 

As  shown  in  Figure  13,  this  system  consisted  of  a  porous, 
water  supply  tube  (positioned  above  the  uppermost  active  tube), 
a  water  heater  which  served  as  a  supply  tank,  a  rotameter  to 


31 


regulate  the  amount  cf  water  supplied  to  the  porous  tube,  a 
pump  driven  by  a  1/2  HP  electric  motor,  and  associated  piping 
and  valves  . 

The  stainless  steel  porous  tube  was  305  mm  long,  the  length 
of  the  test  condenser,  with  non-porous  stainless  steel  tube 
extensions  silver  soldered  on  to  each  end  for  the  purpose  of 
tube  sealing  and  connection  to  the  supply  tubing.  Water  entered 
one  end  of  this  tube;  the  other  end  was  blanked  off. 

The  supply  tank  was  a  25  gallon  hot  water  heater  which  was 
modified  by  the  replacement  of  its  thermostatic  control  with 
a  more  accurate  temperature  controller.  The  temperature  con¬ 
troller  was  capable  of  maintaining  the  water  temperature  to 
within  approximately  1°C  of  its  set  point.  To  aid  in  keeping 
the  active  tabes  free  _ .  contamination  that  could  lead  to 
drop-wise  condensation,  distilled  water  was  used  exclusively 
in  this  system.  The  water  was  pumped  to  the  rotameter  by  the 
1/2  HP  electric  motor  driven  pump,  and  by  regulating  the  rota¬ 
meter  and  the  tank's  recirculation  valve,  the  flow  rate  to  the 
porous  tube  could  be  controlled  with  reasonable  accuracy. 

The  system  was  used  to  supply  water  to  the  top  active  tube 
of  the  test  condenser  at  the  existing  condenser  steam  satura¬ 
tion  temperature.  By  this  method,  the  five  active  tubes  in¬ 
stalled  in  the  test  condenser  could  be  inundated  with  condensate 
from  above  to  simulate  a  tube  column  of  more  than  the  five 
active  tubes . 


32 


V  . 


CONDENSATE  SYSTEM 


The  condensate  system,  as  shown  in  Figure  14 ,  was  composed 
of  the  test  condenser  and  hotwell,  the  secondary  condenser  and 
hotweil,  the  condensate  pump,  and  associated  piping  and  valves. 
The  only  modifications  to  this  sytem  were  those  to  the  test 
condenser  noted  above. 

The  test  condenser  and  secondary  condenser  hotwell  sight 
glasses  were  calibrated  so  that  the  amount  of  condensate  pro¬ 
duced  could  be  measured.  The  condensate  formed  in  the  test 
condenser  was  collected  in  the  test  condenser  hotwell  by  closing 
valve  C-l.  Opening  the  valve  allowed  the  condensate  to  be 
gravity  fed  to  the  secondary  condenser  hotwell.  Additionally, 
condensate  in  the  secondary  hotwell  was  obtained  from  the 
secondary  condenser,  which  condensed  all  steam  not  condensed 
in  the  test  condenser,  and  from  the  desuperheater  drain.  All 
condensate  collected  in  the  secondary  condenser  hotwell  could 
be  removed  by  means  of  an  electrically  driven  condensate  pump 
(1/2  HP)  and  dumped  into  the  building's  drainage  system. 

3y  measuring  the  condensate  produced  per  unit  time  in  the 
test  condenser,  the  porous  tube  rotameter  could  be  set  to 
increment  the  condensate  to  the  top  tube  to  simulate  additional 
tubes.  By  measuring  the  total  amount  of  condensate  produced 
per  unit  time,  in  both  the  test  and  secondary  condensers  and 
substracting  the  amount  of  water  supplied  to  the  porous  tube, 
the  vapor  velocity  could  be  calculated  as  previously  noted. 
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STEAM  SYSTEM 


The  steam  system,  shown  in  Figure  15,  was  unchanged  from 
Morrison's  design.  The  steam  used  for  the  test  apparatus  was 
locally  generated  house  steam.  A  19  mm  O.D.  stainless  steel 
line  supplied  steam  through  a  steam  supply  valve,  MS-3,  to  a 
cost  iron  steam  separator.  Downstream  of  the  separator,  a 
19  mm  O.D.  stainless  steel  line  supplied  steam  to  two  Mupro 
bellows  valves  which  were  used  in  conjunction  with  the  supply 
valve  to  regulate  the  steam  supply  pressure.  From  these  valves, 
the  steam  was  fed  through  another  19  mm  O.D.  stainless  steel 
line  to  the  desuperheater,  and  from  the  desuperheater  through 
a  64  mm  I.D.  stainless  steel  pipe  to  the  test  condenser  diffuser. 
The  steam  supply  pressure  was  monitored  by  a  gage  just  down¬ 
stream  of  the  steam  supply  valve.  A  compound  gage  downstream 
of  the  Mupro  valves  could  be  used  to  monitor  the  pressure  drop 
across  the  valves  if  desired. 

Since  house  steam  was  used,  the  operator  had  no  control 
over  the  state  point,  quality,  non-condensable  gas  content, 
or  contamination  of  the  supply  steam.  Moreover,  it  was  not 
practically  or  economically  feasible  to  monitor  these  condi¬ 
tions  at  the  time  experiments  were  conducted. 

F.  COOLING  WATER  SYSTEM 

The  cooling  water  system  was  a  shown  in  Figure  13.  Normal 
house  water  was  utilized  for  the  test  apparatus,  with  the 
water  stored  in  a  1.2  meter  cubical  plexiglass  tank.  The 
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water  was  pumped  from  the  supply  tank  by  a  5  HP  electrically 
driven  pump  through  a  51  mm  O.D.  plastic  pipe  to  a  header. 

Five  rotameters  were  connected  to  the  header,  with  each  used 
to  control  the  flow  through  a  separate  active  tube.  The  water 
leaving  the  rotameters  passed  through  16  mm  O.D.  stainless 
steel  tubing  to  the  active  tubes.  These  tube  runs  were  at 
least  1.3  m  long  to  ensure  fully  developed  flow  going  into  the 
test  condenser.  The  water  leaving  the  active  tubes  was  returned 
to  the  supply  tank  where  a  7  1/2  HP  electrically  driven  pump  dis¬ 
charged  the  water  through  a  filter  to  a  cooling  tower  to  main¬ 
tain  a  constant  cooling  water  supply  temperature  and  to  keep 
the  water  as  clean  as  possible.  The  cooling  tower  was  located 
outside  the  building  and  consisted  of  four  truck  radiators 
across  which  air  was  blown  by  a  fan.  The  heat  exchanger  and 
fan  were  enclosed  in  a  wooden  structure  with  louvered  openings 
for  ventilation. 

The  tubing  runs  between  the  rotameters  and  active  tubes 
were  divided  into  several  sections.  A  ^5  cm  section  was  con¬ 
nected,  by  5  cm  lengths  of  tygon  tubing  and  hose  clamps,  to 
each  active  tube  and  to  a  12  cm  section  which  housed  the  cool¬ 
ing  water  inlet  thermocouples .  The  12  cm  sections  were  then 
connected,  by  5  cm  lengths  of  tygon  tubing  and  hose  clamps, 
to  the  cooling  water  supply  line.  Thus  on  the  inlet  side,  the 
cooling  water  lines  could  be  readily  disassembled  to  facili¬ 
tate  active  tube  removal  and  installation  or  repairs  to  the 
thermocouples.  On  the  outlet  side  of  the  test  condenser,  12  cm 
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sections  containing  the  outlet  cooling  water  thermocouples 
were  connected,  again  by  tygon  tubing  and  hose  clamps,  to 
the  active  tubes  and  the  discharge  piping. 

G.  VACUUM  SYSTEM 

As  shown  in  Figure  14,  an  air  driven  air  ejector  was 
utilized  to  remove  the  non-condensable  gases  and  to  maintain 
the  test  condenser  vacuum  as  per  Morrison’s  design.  The  air 
ejector  pulled  a  suction  on  the  secondary  condenser  hotwell 
and  discharged  through  a  muffler  to  reduce  the  noise  hazard. 
Air  was  supplied  to  the  air  ejector  at  minimum  of  about 
100  psig  from  an  Ingersol-Rand  air  compressor.  The  system 
was  capable  of  maintaining  a  minimum  test  condenser  vacuum 
of  about  100  mm  Hg  absolute. 

H.  DESUPERHEATER  SYSTEM 

The  desuperheater  system,  shown  in  Figure  15,  was  used  to 
control  the  steam  inlet  temperature  to  the  test  condenser  so 
that  a  maximum  of  10°C  superheat  was  maintained.  This  was 
necessary  to  avoid  a  correction  to  the  heat  transfer  coeffi¬ 
cients  for  h^.  .  In  practice,  the  desuperheater  was  not  re¬ 
quired  for  the  atmospheric  pressure  runs,  but  was  used  for 
the  vacuum  runs  for  which  the  degree  of  superheat  was  con¬ 
trolled  to  about  2°C  of  superheat. 

The  system  consisted  of  a  desuperheater  tank,  a  water 
supply  tank,  a  1/3  HP  electrically  driven  pump,  a  rotameter 
and  associated  piping  and  valves.  The  desuperheater  was 
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constructed  of  a  stainless  steel  tank,  318  mm  in  diameter  and 
503  mm  in  height.  The  tank  top  was  welded  onto  the  discharge 
pipe  and  was  then  bolted  onto  the  tank.  A  rubber  gasket 
provided  for  sealing  between  the  tank  and  the  tank  top.  Four 
fan  type  spray  nozzles  were  inserted  equidistant  around  the 
circumference  and  approximately  10  cm  from  the  top  of  the  de¬ 
superheater  to  supply  the  cooling  water  in  a  fine  spray.  A 
drain  line  was  also  installed  in  the  base  of  the  desuperheater 
tank . 

Cooling  water,  distilled  to  aid  in  active  tube  cleanliness, 
was  stored  in  a  10  gallon  aluminum  tank.  The  water  was  pumped 
from  the  tank  through  9  mm  O.D.  copper  tubing  to  a  rotameter 
and  valve  and  then  to  the  spray  nozzles. 

I.  INSTRUMENTATION 
1 .  Flow  Rates 

a.  The  cooling  water  flow  rate  was  measured  separately 
for  each  active  tube  through  the  use  of  rotameters.  Starting 
from  the  top  active  tube  to  the  bottom  active  tube,  the  supply¬ 
ing  rotameters  had  calibrated  100%  maximum  flow  rates  of  18.5 

t  0.2,  18.5  ±  0.2,  13.5  ±  0.2,  19.2  ±  0.2,  and  17.6  ±  0.2  GPM. 

b.  The  desuperheater  cooling  water  flow  rate  was  mea¬ 
sured  by  a  rotameter  having  a  calibrated  100%  maximum  flow 
rate  of  1242  ±  12  ml/min. 
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c.  The  porous  tube  water  supply  flow  rate  was  measured 
by  a  rotameter  having  a  calibrated  100“  maximum  flow  rate  of 
399  t  10  ml/min. 

d.  All  rotameters  were  calibrated  using  the  procedure 
listed  in  Appendix  A  of  Reference  5  . 

2 .  Temperatures 

Stainless  steel  sheathed  ccpper-constantan  thermo¬ 
couples  were  utilized  as  the  temperature  monitoring  devices  . 
Figure  15  shows  the  locations  of  the  thermocouples.  Two 
thermocouples  measured  the  cooling  water  inlet  temperature, 
and  three  measured  the  cooling  water  outlet  temperature  for 
each  active  tube  in  an  attempt  to  obtain  an  accurate  bulk 
temperature  at  each  sensing  location.  Four  0.51  mm  0.D.  ther¬ 
mocouples  were  silver  soldered  into  axial  grooves  within  the 
wall  of  the  active  tubes  as  noted  above.  Additionally,  thermo¬ 
couples  were  used  to  monitor  the  test  condenser  steam  inlet 
temperature  and  the  temperature  of  the  condensate  in  the  test 
condenser  hotwell.  All  thermocouples  were  calibrated  follow¬ 
ing  the  procedure  of  Appendix  A  of  Reference  5. 

A  Gulton  Industries,  West  23,  0-500°F  temperature  con¬ 
troller  was  used  to  control  the  temperature  of  the  porous  tube 
supply  water.  The  controller  had  a  manufacturers  stated  accu¬ 
racy  of  ±  0.5%  of  the  span  (about  t  1.25°F). 

3 .  Pressure 

a.  The  steam  supply  pressure  was  monitored  by  a 
Bourdon  tube  gage  just  downstream  of  the  steam  supply  valve. 
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b.  The  test  condenser  pressure  was  monitored  by  a 
mercury  manometer  and  by  a  pressure  transducer.  The  pressure 
transducer  was  calibrated  against  a  mercury  manometer. 

c.  Pressure  gages  were  also  used  to  monitor  the  test 
condenser  pressure,  the  secondary  condenser  hotwell  pressure, 
and  the  pressure  downstream  of  the  Mupro  valves  .  These  gages 
were  used  strictly  for  system  operation  purposes. 

4 .  Data  Collection  and  Display 

An  Autodata  Nine  Scanner  was  used  to  record  and  dis¬ 
play  the  thermocouple  and  pressure  transducer  readouts  .  The 
temperatures  were  recorded  in  degrees  Celsius,  while  the  pres¬ 
sure  was  recorded  as  a  voltage  which  was  converted  to  a  reading 
in  mm  Hg  absolute  by  a  calibration  chart.  The  pressure  trans¬ 
ducer  was  assigned  to  channel  1  of  the  Autodata  Nine  Scanner, 
while  the  thermocouples  were  assigned  channels  as  indicated 
in  Table  I. 
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IV  . 


PROCEDURES 


A.  OPERATING  PROCEDURES 

1 .  Tube  Preparation 

Prior  to  installation,  and  periodically  thereafter, 
each  active  tube  was  prepared  in  accordance  with  the 
cleaning  procedure  outlined  in  Appendix  A  of  Reference  8. 

In  addition,  as  a  part  of  the  light-off  procedure  of  the  test 
apparatus,  the  tubes  were  cleaned  by  running  steam  at  atmos¬ 
pheric  pressure  through  the  test  condenser  for  about  twenty 
minutes  without  any  cooling  water  flow  through  the  tubes . 
During  this  steam  cleaning  operation,  the  tubes  were  period¬ 
ically  rinsed  with  distilled  water  utilizing  the  porous  tube 
water  supply  system. 

The  above  mentioned  procedures  were  generally  found 
to  be  sufficient  to  eliminate  any  visible  signs  of  drop-wise 
condensation  on  the  active  tubes.  When  drop-wise  condensation 
did  occur  during  the  progress  of  a  day's  set  of  runs,  it  was 
discovered  that  simply  rinsing  the  tubes  using  the  porous  tube 
water  supply  system  was  sufficient  to  restore  film-wise 
condensation . 

2 .  System  Operation  and  Steady  State  Conditions 

Complete  operating  instructions  are  listed  in  Appendix 

A  and  Appendix  B . 


40 


A  steady  state  condition  was  considered  achieved  when 
the  cooling  water  inlet  temperature  did  not  vary  more  than 
±  0 ,5°C/HR  and  the  steam  inlet  temperature  did  not  vary  more 
than  ±  0.3°C/min.  In  general,  a  steady  state  condition  was 
reached  about  two  hours  after  initial  system  light-off  and 
about  fifteen  to  thirty  minutes  after  changes  to  the  cooling 
water  or  porous  tube  supply  water  flow  rates  . 

Once  a  steady  state  condition  was  obtained,  the  runs 
were  made.  The  length  of  each  run,  or  each  increment  of  five 
tubes,  was  five  minutes.  The  following  data  were  taken  at  one 
minute  intervals  for  each  run: 

a.  the  test  condenser  and  secondary  condenser 
hotwell  levels, 

b.  the  manometer  reading, 

c.  the  steam  supply  pressure, 

d.  the  setting  of  each  rotameter, 

e.  the  thermocouple  readings,  and 

f.  the  pressure  transducer  reading 

B.  DATA  REDUCTION  PROCEDURES 

To  simplify  the  data  reduction,  actual  tube  wall  tempera¬ 
tures  were  obtained.  Utilizing  the  wall  temperatures  directly 
eliminated  the  necessity  of  using  the  Wilson  Plot  technique 
for  calculating  the  steam  side  heat  transfer  coefficient. 

The  following  approach,  that  of  using  standard  heat  transfer 
equations  and  an  energy  balance,  was  used  to  evaluate  the 
raw  data: 
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1.  Total  Heat  Transferred  Der  Tube 


=  m- c  •  ( T  -  T  . ) 

11  CO  Cl 


(10) 


where , 


m. 

i 


(T  -T  ■ ) . 
co  ci  i. 


=  total  heat  transferred  for  tube  i  per  unit 
time  CkW] 

=  cooling  water  mass  flow  rate  for  tube  i 
[kg/s ] 

=  the  difference  between  the  outlet  and  inlet 
cooling  water  temperature  for  tube  i  [°C] 
c^  =  the  specific  heat  of  the  cooling  water 
[kJ/kg  *  °Cl 

Steam  Side  Heat  Transfer  Coefficient 


hi 


=  <U/Ai  (Ts  -  V 


(11) 


where  , 


hi 


Ai 


T 

w 


=  the  steam  side  heat  transfer  coefficient  for 
tube  i  [kW/m^  •  °C] 

2 

=  the  heat  transfer  area,  m 
=  the  steam  saturation  temperature  [°C] 

=  the  average  tube  wall  temperature  [°C] 

Sample  calculations  are  presented  in  Appendix  D. 

The  steam  side  heat  transfer  coefficient  was  the  para¬ 
meter  used  to  determine  the  performance  of  a  tube  in  the  tube 
bundle.  Three  assumptions  were  made  in  applying  Equation  (11) 
to  the  data  obtained  from  the  experiment: 
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1.  The  resistance  due  to  non-condensable  gases  was  negli¬ 
gible  .  The  system  was  checked  for  tightness  and  was 
found  to  be  satisfactory.  Mon-condensable  gases  brought 
in  by  the  steam  were  exhausted  by  the  air  ejector  as 
the  good  steam  flow  conditions  in  the  test  condenser 
continually  swept  non-condensable  gases  away  from  the 
active  tubes  . 

2.  The  resistance  of  the  copper  active  tubes  was  considered 
negligible.  The  calculated  1°C  AT  across  the  tube 
wall  was  accounted  for  by  assigning  a  significant  un¬ 
certainty  to  the  average  tube  wall  temperature  which 
was  determined  from  the  thermocouples. 

3.  Subcooling  of  the  condensate  film  on  the  tube  walls 
was  neglected. 

C.  DATA  REDUCTION  PROGRAM 

A  computer  program  was  utilized  to  analyze  the  data.  The 
program  was  in  Basic  language  and  was  run  on  an  HP- 35  computer 
system.  A  peripheral  plotter  was  used  to  plot  the  results. 

The  program  is  presented  in  Appendix  E. 


V.  RESULTS  AND  DISCUSSION 


A.  3ACKGR0UND  INFORMATION 

Sixty-two  runs  were  made.  Runs  1-40  and  Run  62  were  con¬ 
ducted  with  the  test  condenser  at  atmospheric  pressure .  Steam 
supply  pressures  of  50  psig  (Runs  1-10  and  Run  62),  45  psig 
(Runs  11-20),  40  psig  (Runs  21-30),  and  35  psig  (Runs  31-40) 
were  utilized.  To  ensure  that  non-condensable  gases  were 
flushed  out  of  the  test  condenser,  the  above  runs  were  con¬ 
ducted  with  the  steam  exhaust  pipe  from  the  test  condenser  to 
the  secondary  condenser  disconnected  and  replaced  by  a  1.5  cm 
O.D.  pipe  which  dumped  the  exhaust  steam  into  the  bilges  (as 
shown  in  Figure  7a).  Runs  41-61  were  conducted  under  vacuum 
at  steam  supply  pressures  of  35  psig  (Runs  41-45),  40  psig 
(Runs  46-50),  45  psig  (Runs  51-55),  and  50  psig  (Runs  56-61). 
All  runs  were  made  with  a  cooling  water  mass  flow  rate  of 
10.98  kg/min  for  each  active  tube.  The  desuperheater  was 
required  for  only  the  vacuum  runs,  for  which  a  cooling  water 
volumetric  flow  rate  of  100  ml /min  was  used.  Runs  1-60  were 
made  with  no  water  being  supplied  through  the  porous  tube. 
Using  the  porous  tube  water  supply  system  to  simulate  addi¬ 
tional  tubes,  Runs  61  and  62  were  conducted  for  tube  bundles 
of  30  and  10  tubes  deep,  respectively. 

The  data  collected  for  each  run  are  shown  in  Tables  II 
through  IV.  It  should  be  noted  that  the  data  shown  are  Che 
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averages  of  the  five  data  sets  taken  for  each  run,  corrected 
using  thermocouple,  pressure  transducer,  rotameter,  and  hot- 
well  volume  calibration  curves.  A  full  set  of  data  was 
collected  for  the  vacuum  runs  ,  but  due  to  the  dumping  of 
exhaust  steam  into  the  bilges  during  the  atmospheric  pressure 
runs ,  the  volumetric  flow  rate  of  condensate  from  the  second¬ 
ary  condenser  could  not  be  measured.  This  also  prohibited 
vapor  velocity  calculation. 

The  calculated  results  are  tabulated  in  Table  V.  The 
results  are  presented  to  the  appropriate  number  of  significant 
dig.its  as  determined  from  the  uncertainty  analysis  presented 
in  Appendix  D. 

B.  OBSERVATIONS 

Film-wise  condensation  without  any  visible  evidence  of 
drop-wise  condensation  was  observed  for  all  tubes  during  all 
runs.  This  was  accomplished  by  the  cleaning  procedures  de¬ 
tailed  earlier.  It  must  be  noted  that  only  the  front  half  of 
each  tube  could  be  viewed  through  the  test  condenser  window, 
therefore  drop-wise  condensation  could  have  occurred  unnoticed 
on  the  back  halves  of  the  tubes .  Based  on  the  observation  of 
the  active  tubes  and  the  tube  wall  temperature  data,  the 
effect  of  drop-wise  condensation  was  considered  to  be  negligible. 

Observation  of  the  motion  of  the  condensate  drops  during 
both  atmospheric  pressure  and  vacuum  conditions  revealed  a 
general  logitudinal  drop  migration  on  all  tubes.  The  migration 


was  left  to  right  under  atmospheric  pressure  conditions  and 
right  to  left  under  vacuum  conditions.  During  several  opera¬ 
tions  requiring  the  disassembly  of  the  test  condenser  flanges, 
it  was  noted  that  the  piping  to  and  from  the  condenser  did  not 
align  properly  such  that  piping  had  to  be  forced  into  position 
for  flange  alignment  and  assembly.  Apparently,  when  the  steam 
exhaust  pipe  was  disconnected  for  atmospheric  pressure  opera¬ 
tion,  the  orientation  of  the  test  condenser  changed, resulting 
in  a  change  in  the  direction  of  drop  migration.  By  adjustment 
of  the  leveling  nuts  on  the  test  condenser  support  bracket, 
the  amount  of  drop  migration  was  reduced  to  a  level  at  which 
its  affect  was  considered  negligible.  It  was  also  noted  that 
at  certain  tube  localities,  drops  moved  randomly  either  left 
or  right  possibly  due  to  localized  tube  irregularities  since 
this  observation  was  not  widespread.  All  evidences  of  drop 
migration  were  then  believed  to  have  resulted  from  condenser 
leveling  and  tube  irregularity  problems ,  and  not  the  result  of 
secondary  vapor  flow. 

The  spacing  of  the  tubes,  24  mm  from  tube  center  to  tube 
center  or  8  mm  between  tubes,  led  to  several  interesting  con¬ 
densate  drop  observations.  Under  of  influence  of  low  conden¬ 
sate  production  or  inundation  rates ,  the  drops  from  an  upper 
tube  would  break  away  and  fall  less  than  two  millimeters 
before  striking  the  tube  below.  As  a  consequence,  no  splash¬ 
ing  occurred  and  the  drops  appeared  to  spread  in  all  directions 
radially  from  the  point  of  drop  impact,  incorporating  into 
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the  condensate  film  rather  than  rolling  around  the  tube  peri¬ 
meter.  At  higher  condensate  rates,  such  as  under  atmospheric 
pressure  or  under  inundation  from  the  porous  tube  for  vacuum 
runs ,  the  condensate  drops  from  the  upper  tube  elongated  all 
the  distance  to  the  next  lower  tube.  The  drops  instantaneously 
joined  the  upper  and  lower  tube  and  had  the  appearance  of  being 
sucked  into  the  condensate  film  on  the  lower  tube.  It  was  also 
noted  that  as  the  condensation  or  inundation  rate  increased, 
drops  formed  at  more  locations  along  a  given  tube  but  that  the 
ratios  of  drops  per  tube  between  the  tubes  remained  nearly 
constant . 

It  is  interesting  to  note  that  Shklover  and  3uevich  [Ref. 
17]  conducted  studies  of  drop  behavior  for  condenser  tubes 
spaced  3  mm  apart.  Though  their  operating  conditions  were 
not  the  same  as  used  in  this  study,  they  observed  similar 
behaviors  as  described  above,  except  that  they  observed  the 
drops  to  roll  over  the  tube  surface  and  did  not  observe  the 
case  in  which  the  drop  separated  from  the  upper  tube  before 
reaching  the  lower  tube . 

C.  AVERAGE  TUBE  WALL  TEMPERATURES 

As  noted  previously,  the  use  of  Equations  (10)  and  (11) 
to  calculate  the  steam  side  heat  transfer  coefficient  was 
made  possible  by  the  installation  of  thermocouples  90°  apart 
within  the  tube  walls.  However,  the  questions  arose  as  to 
how  accurately  the  average  tube  wall  temperature  could  be 


measured  and  how  the  average  tube  wall  temperature  was  to 
be  defined  for  the  direct  wall  temperature  measurement. 

Figure  17  is  a  polar  plot  of  a  typical  tube  wall  tempera¬ 
ture  distribution.  It  can  be  seen  that  the  distribution  is 
asymmetric  with  the  majority  of  the  tube  surface  area  at  a 
temperature  higher  than  the  average  of  the  top  and  bottom  wall 
temperatures.  This  is  a  consequence  of  the  nonuniform  nature 
of  the  condensate  film  thickness  around  the  tube,  with  the 
film  thickening  from  top  to  bottom.  The  film  tends  to  insulate 
the  wall,  yielding  lowest  wall  temperatures  where  the  film  is 
thickest . 

The  temperature  distribution  becomes  even  more  complex  when 
events  which  locally  thicken  or  thin  the  condensate  film  are 
considered.  A  drop  forming  at  the  base  of  a  tube  or  falling 
onto  the  top  of  the  tube  will  cause  temporary  reduction  in  the 
wall  temperature  at  the  point  of  occurrence.  A  small  area  of 
drop-wise  condensation  or  a  localized  turbulent  eddy  can  locally 
increase  the  wall  temperature.  When  side  drainage  occurs,  a 
tube  can  have  its  condensate  film  either  thinned  or  thickened 
resulting  in  a  higher  or  lower  wall  temperature,  respectively. 
Additionally,  high  vapor  velocities  can  either  strip  away  part 
of  the  film  layer  or  cause  increased  drop  splashing  or  rippling 
of  the  film  all  of  which  increase  the  tube  wall  temperature. 

The  above  description  demonstrates  that  the  best  approxima¬ 
tion  of  the  average  wall  temperature  would  be  attained  when 
the  temperature  was  measured  at  as  many  points  as  possible  over 
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the  entire  tube  wall  and  for  an  extended  period  of  time.  Then, 
an  average  of  all  the  data  would  be  needed.  The  tube  size, 
method  of  thermocouple  installation,  and  the  cost  of  the 
thermocouples  limited  the  number  of  points  of  temperature  mea¬ 
surement  to  four.  Since  the  tube  wall  temperature  would  also 
vary  through  the  wall  thickness,  it  was  decided  to  install  the 
four  thermocouples  90  degrees  apart  (top,  bottom,  and  two  sides) 
in  grooves,  as  noted  previously,  such  that  the  thermocouple 
sensing  point  would  be  at  the  midpoint  of  the  wall  thickness. 
Calculations  made  using  the  data  obtained  from  the  test  runs 
revealed  that  the  AT  across  the  tube  wall  was,  on  the  average, 
about  1°C.  Therefore,  fairly  small  deviations  in  thermocouple 
location  could  also  have  a  significant  effect  on  the  calculated 
average  tube  wall  temperature . 

Considering  all  the  above,  and  to  avoid  giving  undue  weight 
to  the  two  side  wall  temperatures,  it  was  decided  that  the  best 
average  tube  wall  temperature  would  be  obtained  by  averaging 
the  side  wall  temperatures,  then  averaging  this  result  with 
the  top  and  bottom  wall  temperatures.  It  was  also  decided 
that  the  average  tube  wall  temperature  must  be  considered  to 
be  subject  to  fairly  significant  uncertainty.  Additionally, 
due  to  damage  to  thermocouples  and  the  inability  to  replace 
them  at  the  time  tests  were  conducted,  only  the  center  active 
tube  had  a  full  compliment  of  thermocouples.  Tubes  1,  4,  and 
5  had  three  thermocouples  each  (top,  bottom,  and  side),  while 
tube  2  had  only  two  (top  and  bottom).  For  these  partially 
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ir-strumented  tubes ,  the  average  tube  wall  temperature  was  cal¬ 
culated  by  averaging  the  readings  from  the  available  thermo¬ 
couples.  It  should  be  noted  that  this  situation  was  most 
disadvantageous  for  the  case  of  tube  2  since  the  wall  tempera¬ 
ture  obtained  for  this  tube  was  certainly  lower  than  would  be 
measured  with  the  side  thermocouples  available.  In  this 
regard,  it  can  be  noted  from  the  tube  wall  temperatures  tab¬ 
ulated  in  Table  II,  that  the  second  tube  consistently  had  an 
average  wall  temperature  much  below  that  expected.  Thus,  to 
be  conservative,  an  uncertainty  of  ±  1.0°C  was  assigned  to  all 
wall  temperatures . 

We  now  turn  to  two  other  important  considerations  in  regard 
to  the  average  wall  temperatures.  These  are  the  effect  of 
vapor  velocity  on  the  average  tube  wall  temperature  and  the 
consideration  of  how  the  location  of  a  tube  in  a  tube  bundle 
affects  its  average  wall  temperature. 

The  average  tube  wall  temperature  should  increase  as  the 
vapor  velocity  increases .  This  is  the  result  of  increased 
vapor  shear  which  can  thin  the  condensate  film  by  causing 
small  portions  of  the  film  to  be  removed  from  the  tube  and 
entrained  in  the  vapor  and  by  increasing  the  velocity  of  the 
condensate  flowing  on  the  tube  .  For  the  vacuum  condition 
runs,  where  the  mean  vapor  velocity  was  calculated,  the  re¬ 
sults  shown  in  Table  VI  confirm  the  increase  in  the  average 
tube  wall  temperature  with  increased  mean  vapor  velocity  even 
for  the  relatively  small  vapor  velocity  increases  indicated. 
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This  could  be  seen  for  all  five  tubes  of  the  tube  bundle,  but 
was  most  apparent  for  the  first  two  tubes.  Mote  also  that  the 
increase  in  the  average  tube  wall  temperature  resulted  in  a 
dropping  off  of  T  -  T  which  was  manifested  in  a  higher  heat 
transfer  coefficient  with  increased  mean  vapor  velocity. 

Figures  18  and  19  depict  the  decrease  in  the  average  tube 
wall  temperature  through  a  tube  bundle.  For  the  vacuum  condi¬ 
tion  run  (Run  61)  and  the  atmospheric  pre^  ure  run  (Run  62), 
the  general  trend  was  one  of  a  sharp  initial  drop  in  the  average 
tube  wall  temperature  over  the  upper  tubes  followed  by  a  much 
more  gradual  reduction  through  the  remainder  of  the  tube  bundle. 
This  was  best  illustrated  for  Run  51  where  about  70%  of  the 
total  wall  temperature  drop  occurred  in  the  first  ten  tubes . 
These  drops  in  the  average  tube  wall  temperature  through  a 
tube  bundle  can  be  attributed  to  a  decreasing  vapor  velocity 
as  the  mass  of  the  steam  is  reduced  by  condensation  through 
the  tube  bundle  and  to  an  increasing  amount  of  condensate  inun¬ 
dation.  Both  of  these  factors  would  serve  to  lower  the  average 
tube  wall  temperatures  by  allowing  for  an  increase  in  the  con¬ 
densate  film  thickness.  It  must  be  considered,  however,  that 
the  above  results  were  obtained  through  the  use  of  the  porous 
tube  water  supply  system,  the  operation  of  which  still  requires 
refinement . 

A  study  of  Table  II  as  well  as  Figures  18  and  19  demon¬ 
strated  that  the  average  tube  wall  temperatures  for  the  second 
and  fifth  tubes  of  the  test  condenser  were  lower  than  would 
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be  exoected. 


This  discrepancy  for  the  second  tube  was  due  at 


least  in  part  to  the  lack  of  functioning  thermocouples  on  both 
sides  of  the  tube.  Additional  test  runs  with  complete  instru¬ 
mentation  on  this  tube  will  be  required  before  further  extlan- 
ation  can  be  provided.  In  regard  to  the  fifth  tube,  it  must 
be  noted  that  the  vapor  velocity  across  this  tube  will  be  re¬ 
duced  by  the  effect  of  an  increase  in  the  steam  flow  area  below 
the  tube  due  to  the  absence  of  additional  tubes.  This  could 
significantly  reduce  the  tube  wall  temperature.  It  is  there¬ 
fore  considered  advisable  to  install  a  row  of  dummy  tubes 
below  the  fifth  tube  for  test  runs  in  which  additional  tubes 
are  to  be  simulated.  However,  these  dummy  tube  should  be 
removed  for  runs  in  which  additional  tubes  are  not  to  be 
simulated . 

D.  SINGLE  TUBE  RESULTS 

One  method  of  gauging  the  validity  of  experimental  results 
is  to  make  comparisons  with  theoretical  predictions  and  with 
experimental  correlations.  If  one  compares  the  results  for 
only  the  top  tube  of  the  five  active  tube  column  used  for  the 
test  runs  with  such  predictions  and  correlations,  then  one 
cannot  only  determine  the  validity  of  the  top  (single)  tube 
results,  but  also  obtain  a  better  insight  into  the  validity 
of  the  multiple  tube  results. 

Tables  VII  and  VIII  compare  the  experimentally  obtained 
heat  transfer  coefficient  (h^)  to  the  Nusselt  prediction 
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(h„  )  and  to  the  Fujii  correlation  (hQ)  .  Remember  that  h^ 
was  calculated  from  Equation  (11),  h , . u  from  Equation  (1), 
and  h  by  using  Equations  (9)  and  (12). 

Table  VII  displays  the  results  for  the  atmospheric  runs, 
with  hQ  not  calculated  since  vapor  velocities  could  not  be 
determined.  It  should  be  noted  that  although  the  vapor  velo¬ 
cities  were  unavailable,  the  vapor  velocities  should  be  de¬ 
creasing  from  a  highest  value  for  Runs  1-10  to  a  lowest  value 
for  Runs  31-40  in  response  to  decreasing  steam  supply  pressure 
and  therefore  decreasing  system  pressure  drops.  Observation 
of  the  exhaust  steam  dumped  into  the  bilges  confirmed  the 
above  and  also  confirmed  the  fact  that  fairly  quiescent  steam 
conditions  prevailed. 

The  results  showed  that  the  experimental  and  theoretical 
values  for  the  heat  transfer  coefficients  were  compatible, 
with  the  experimental  values  on  the  order  of  10%  higher  than 
the  Musselt  theory  prediction.  This  was  within  the  as  much  as 
20%  discrepancy  obtained  in  many  studies,  and  could  be  accounted 
for  by  the  difference  between  the  conditions  assumed  by  the 
Musselt  theory  and  the  actual  conditions  of  condensation  on 
a  horizontal  tube.  Moreover,  it  must  be  considered  that  the 
single  tube  under  study  was  not  a  single  tube  at  all  since 
the  tube  was  positioned  between  the  porous  tube  and  the  second 
active  tube.  Knudsen  and  Katz  [Ref.  22]  ha;e  shown  that  for 
air  flowing  downward  over  a  staggered  tube  bundle,  the  first 
tube  can  cause  turbulence  which  raises  the  local  Musselt 
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number  (and  therefore  the  local  heat  transfer  coefficient) 
for  the  tubes  beneath  it  with  the  effect  being  greatest  for 
the  second  and  ruird  tubes.  Thus,  the  porous  tube  could 
cause  an  increase  in  the  heat  transfer  coefficient  for  the 
top  active  tube.  Additionally,  when,  as  observed,  condensate 
drops  do  not  separate  from  an  upper  tube  but  rather  form  a 
link  between  an  upper  and  lower  tube,  it  may  be  that  more 
condensate  is  removed  from  the  upper  tube  than  otherwise  would 
take  place.  The  result  may  be  that  the  tube  below  the  studied 
"single  tube"  assisted  the  "single  tube"  and  raised  its  heat 
transfer  coefficient,  though  this  is  an  area  which  requires 
further  study  . 

The  results  also  indicated  the  lack  of  a  clear,  distinctive 
relationship  between  h-^ ,  h^u,  and  the  mean  vapor  velocity  trend 
for  the  atmospheric  runs.  Apparently,  the  degree  of  change  in 
the  mean  vapor  velocity  over  the  runs  was  not  sufficient,  under 
the  noted  condition  of  high  inundation,  to  cause  a  marked  effect 
on  the  heat  transfer  coefficients. 

Table  VIII  presents  the  results  for  the  vacuum  runs.  The 
experimental  heat  transfer  coefficient  values  (h-^)  are  seen 
to  be  considerably  higher  than  predicted  by  flusselt  (hNu} ’ 
with  the  extreme  case  being  Run  58  where  h^  was  75%  higher 
than  the  corresponding  h^u.  Upon  comparing  h-^  to  hQ ,  the 
later  of  which  considers  the  effect  of  vapor  velocity,  the 
results  were  seen  to  closely  agree,  especially  for  the  runs 
at  higher  vapor  velocities.  Thus,  it  was  apparent  that  the 
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vapor  velocity  effects  were  significant  for  the  vacuum  runs. 
Mote  also  that  as  the  vapor  velocity  was  increased,  h^u  in¬ 
creased  only  moderately,  while  h-,  and  hQ  increased  dramatic¬ 
ally  which  indicated  that  vapor  shear  was  more  effective  in 
increasing  the  heat  transfer  coefficient  as  the  vapor  velocity 
increased . 

A  comparison  of  the  heat  transfer  coefficients  for  the 
atmospheric  pressure  and  vacuum  runs  revealed  that  both  h^ 
and  h^  values  were  consistently  greater  for  the  vacuum  runs. 

It  was  also  noted  that  the  difference  in  the  h^  values  between 
the  two  sets  of  runs  was  much  greater  than  the  difference  in 

the  h.T  values . 

Nu 

There  are  two  major  factors  which  accounted  for  the  above 
comparisons.  First,  as  has  been  noted,  the  effect  of  the  vapor 
velocity  was  greater  for  the  vacuum  runs .  Increased  vapor 
velocity  increases  the  vapor  shear,  causing  a  thinning  of  the 
condensate  film  on  the  tube,  which  results  in  higher  average 
tube  wall  temperatures.  Noting  that  both  the  equations,  the 
equations  for  h-,  and  h,t  ,  have  a  CT  -  T  )  term  in  their  denom- 
inators ,  raising  the  tube  wall  temperature  raises  the  heat 
transfer  coefficient,  though  to  a  much  higher  degree  for  the 
h^  calculation.  Secondly,  the  rate  of  condensation,  and 
therefore  the  thickness  of  the  tube  wall  condensate  film, 
was  much  greater  for  the  runs  at  atmospheric  pressure.  This 
can  be  seen  by  comparing  the  condensation  rates,  ,  in 
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Table  III.  For  the  atmospheric  pressure  runs,  condensate  pro¬ 
duction  for  the  complete  tube  bundle  averaged  530  ml/min  com¬ 
pared  to  only  260  ml/min  for  the  vacuum  runs.  The  greater 
condensate  film  thickness  for  the  atmospheric  runs  resulted 
in  lowered  tube  wall  temperatures  and  therefore  lowered  heat 
transfer  coefficients,  relative  to  the  vacuum  runs.  This  could 
be  seen  by  comparing  the  relative  differences  between  the  sat¬ 
uration  temperatures  and  average  tube  wall  temperatures  for 
the  atmospheric  pressure  runs  (on  the  order  of  a  25%  difference) 
and  for  the  vacuum  runs  (on  the  order  of  a  15%  difference)  in 
Table  II.  Noting  again  that  the  calculation  of  h-^  is  more 
sensitive  to  a  change  in  (T  -  Tw) ,  the  h^  values  increased 
much  more  between  atmospheric  pressure  and  vacuum  runs  than 
did  the  hjj  values . 

In  summary,  the  heat  transfer  coefficients  for  the  top 
active  tube  were  in  close  agreement  to  the  Nusselt  theory  for 
the  atmospheric  runs,  but  were  far  above  the  Nusselt  theory 
and  in  close  agreement  with  the  Fujii  correlation  for  the 
vacuum  runs .  This  can  be  explained  by  the  dominance  of  the 
effect  of  a  thick  condensate  film  for  the  atmospheric  runs 
and  by  the  dominance  of  vapor  velocity  effects  for  the 
vacuum  runs . 

E.  MULTIPLE  TUBE  RESULTS 

Table  V  and  Figures  20-35  display  the  results  for  the 
multiple  tube  runs.  Due  to  the  large  number  of  runs,  the 
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normalized  average  and  local  heat  transfer  coefficients, 
n Xj / h ^  and  h^/h-^,  respectively,  were  averaged  for  runs  at  the 
same  steam  supply  and  test  condenser  pressures  giving  averaged 
normalized  coefficients,  (E^/h^)  and  (hXj/h^)  .  These  were  then 
plotted  versus  the  tube  number  for  presentation  in  the  Figures. 
The  averaging  of  results  was  considered  justified  in  view  of 
the  close  agreement  of  the  results  at  the  same  pressures  and 
in  view  of  the  reduction  in  the  uncertainty  of  the  results 
afforded  by  the  averaging  process.  The  uncertainty  analysis 
is  given  in  Appendix  C. 

The  data  points  plotted  in  the  Figures  were  fit,  by  the 
least  squares  method,  to  the  Nusselt  type  equations: 


and , 


For  each  set  of  runs,  the  following  values  for  S  and  P  were 
obtained : 


RUN  # 

S 

P 

1-10 

0.21 

0  .78 

11-20 

0.19 

0  .  80 

21-30 

0  .22 

0.78 

31-40 

0.21 

0  .79 

41-45 

0  .12 

0  .85 

46-50 

0.14 

0  .85 

51-55 

0.13 

0.86 

56-60 

0  .17 

0.32 
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Mote  that  for  the  Musselt  equations, 

3  =  0.25  and  P  =  0.75. 

A  general  trend  seen  in  the  Table  and  Figures,  for  both 
atmospheric  pressure  and  vacuum  runs,  was  that  the  heat  trans¬ 
fer  coefficient,  the  normalized  heat  transfer  coefficients, 
and  the  heat  flux  all  decreased  through  the  tube  bundle  for 
each  run  or  set  of  runs.  This  demonstrated  the  effect  of  con¬ 
densate  inundation  through  the  tube  bundle.  The  largest  decrease 
in  both  h  and  h^/h-^  relative  to  the  previous  tube  occurred  for 
the  second  tube  with  the  next  largest  decrease  for  the  fifth 
tube,  while  decreases  for  the  third  and  fourth  tubes  were 
notably  less.  The  h^/h-^  values  showed  the  largest  decrease 
relative  to  the  previous  tube  for  the  second  tube  followed  by 
a  gradual  decrease  through  the  remaining  tubes.  Mote  that  the 
h^/hi  plot  smoothed  out  some  of  the  irregularity  seen  in  the 
correponding  h^/h-^  plot  for  the  same  run.  Also  for  both 
atmospheric  pressure  and  vacuum  runs,  the  heat  transfer  coef¬ 
ficient  for  a  given  tube  increased  as  the  vapor  velocity  in¬ 
creased,  with  the  effect  being  more  pronounced  for  the  vacuum 
runs  . 

Interestingly,  the  normalized  local  and  average  heat  trans¬ 
fer  coefficients  for  a  given  tube  decreased  with  increased 
vapor  velocity  for  the  vacuum  runs.  No  similar  trend  existed 
for  the  atmospheric  pressure  runs.  This  was  most  clear  when 
the  (h^j/h^)  and  (h^/h^)  values  were  compared  for  a  given  tube 
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in  Table  V.  These  results  are  contrary  to  the  expected  trend 
in  that  they  indicate  that  as  the  vapor  velocity  increases, 
the  experimentally  determined  normalized  heat  transfer  coef¬ 
ficients  more  closely  agree  with  Nusselt's  theory.  This  can 
be  seen  from  the  Figures  and  from  ti.e  values  for  the  S  and  P 
exponents  tabulated  above. 

Actually,  the  relationship  between  vapor  velocity  and  the 
normalized  heat  transfer  coefficients  given  in  the  previous 
paragraph  is  extremely  misleading.  One  is  led  to  believe  that 
an  increase  in  vapor  velocity  degrades  the  overall  performance 
of  the  condenser,  when  in  fact  the  heat  flux  (q/A)  and  the 
heat  transfer  coefficient  (h)  for  both  a  given  tube  and  the 
tube  bundle  as  a  whole  increase  as  the  vapor  velocity  increases. 

The  inconsistency  in  the  results  is  resolved  when  one  considers 
the  fact  that,  as  shown  in  Table  V,  the  heat  transfer  coeffi¬ 
cient  for  the  first  tube  (h-^)  increased  notably  more  with  an 
increase  in  vapor  velocity  than  did  the  heat  transfer  coeffi¬ 
cients  for  the  lower  four  tubes  of  the  bundle.  (This  could  be 
due  to  turbulent  steam  flow  over  the  first  tube  caused  by  the 
porous  tube  above  it  and/or  due  to  a  decrease  in  the  vapor 
velocity  through  the  test  condenser  as  the  steam  is  condensed.) 

Then  when  and  especially  h^  are  normalized  by  dividing  by 
h^,  the  resulting  hpj/h^  and  h^/h^  values  naturally  decreased 
as  the  vapor  velocity  increased.  It  is  therefore  apparent 
that  though  the  normalized  heat  transfer  coefficients  are 

j 


59 


invaluable  in  studying  the  effects  of  condensate  inundation 
through  a  tube  bundle,  they  can  cause  confusion  and  faulty 
interpretation  of  results  when  used  to  consider  the  effects 
of  vapor  velocity. 

One  last  comparison  between  the  atmospheric  pressure  and 
vacuum  runs  needs  to  be  made .  It  was  determined  that  smaller 
relative  decreases  in  h,  fT^/h^,  and  h^/h^  occurred  for  each 
tube  through  the  tube  bundle  for  vacuum  condition  runs  than 
for  the  atmospheric  pressure  runs.  This  result  was  attributed 
to  the  greater  inundation  of  tubes  for  the  atmospheric  runs 
and  the  greater  vapor  velocities  for  the  vacuum  runs.  In  other 
words,  the  effects  of  vapor  velocity  to  increase  the  above  men¬ 
tioned  coefficients  dominated  over  the  effects  of  condensate 
inundation  to  lower  them  for  the  vacuum  runs,  while  just  the 
reverse  was  true  for  the  atmospheric  pressure  runs . 

The  final  area  of  multiple  tube  study  was  to  investigate 
the  use  of  the  porous  tube  water  supply  system  to  simulate  a 
tube  bundle  with  more  than  five  active  tubes  in  a  column.  A 
major  concern  was  that  the  average  tube  wall  temperatures 
should  fail  from  the  upper  to  lower  tubes  as  would  be  expected 
for  an  actual  tube  bundle  with  N  number  of  tubes  in  a  column. 
Table  II  and  Figure  18  revealed  that  this  condition  was  satis¬ 
fied  for  the  first  ten  tubes  of  Run  61,  but  that  thereafter, 
the  top  tube  of  each  five  tube  increment  was  too  high.  Table 
II  and  Figure  19  demonstrated  that  for  Run  62  the  sixth  tube 


had  a  higher  tube  wall  temperature  than  the  fifth  tube.  The 
effect  of  the  average  tube  wall  temperature  discrepancies  was 
reflected  in  the  heat  transfer  coefficients,  as  shown  in  Table 
V  and  Figures  36-39,  which  did  not  fall  as  expected,  but  rather 
followed  the  trend  of  the  average  tube  wall  temperatures.  The 
wall  temperature  discrepancies  were  not  excessive;  therefore, 
it  was  believed  that  a  correction  could  be  accomplished  by 
setting  the  porous  tube  inundation  water  temperature  to  the 
condensate  temperature  as  measured  coming  off  the  bottom  tube 
of  the  previous  five  tube  increment.  This  would  allow  for  the 
effect  of  condensate  subcooling  which  was  apparently  more  signi¬ 
ficant  than  originally  anticipated. 

Taking  the  above  into  account,  a  study  of  Figures  36-39  can 
give  a  good  indication  of  the  effect  of  condensate  inundation 
for  larger  depth  tube  bundles.  Run  61,  a  vacuum  run,  demon¬ 
strated  that  a  least  square  curve  fit  to  a  larger  number  of 
tubes  was  much  closer  to  the  Nusselt  prediction  than  for  the 
case  of  a  five  tubes  per  column  bundle.  It  also  appeared 
that  the  vapor  velocity  effects  were  dominant  over  the  first 
four  or  five  tubes,  while  the  effects  of  condensate  inundation 
prevailed  for  the  remainder  of  the  tube  bundle.  Run  62,  an 
atmospheric  pressure  run,  did  not  show  an  appreciable  change 
in  heat  transfer  coefficient  relationships  with  the  addition 
of  tubes  as  compared  to  the  five  tubes  per  column  runs. 

The  least  squares  curve  fits  for  Runs  61  and  62  yielded 
exponents  as  follows : 
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RUN  * 


o 


c 

61  0.23  0.75 

62  0.21  0.78 

From  the  Figures,  it  was  noted  that  the  curve  fits  were  poor 
over  the  upper  portion  of  the  tube  bundles,  but  were  good  for 
the  lower  portion. 
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VI.  CONCLUSIONS 


As  a  result  of  the  above  mentioned  tests,  the  following 
conclusions  were  reached: 

1.  The  average  tube  wall  temperatures  increased  as  the 
vapor  velocity  increased,  but  decreased  as  the  amount 
of  condensate  inundation  increased.  The  heat  transfer 
coefficients  followed  this  same  trend. 

2.  The  experimentally  obtained  values  for  the  heat  trans¬ 
fer  coefficient  and  the  normalized  average  and  local 
heat  transfer  coefficients  were  compatible  with  the 
Nusselt  predictions  for  the  runs  conducted  at  atmos¬ 
pheric  pressure.  The.  same  coefficients  for  the  runs 
conducted  under  vacuum  conditions  greatly  exceeded 
the  Nusselt  predictions  due  largely  to  the  effects 

of  increased  vapor  velocity.  For  the  vacuum  runs, 
the  heat  transfer  coefficient  for  the  top  tube  of  the 
tube  bundle  agreed  closely  with  the  experimental  cor¬ 
relation  proposed  by  Fujii. 

3.  The  simulation  of  a  tube  bundle  of  more  than  five 
active  tubes  per  column  under  vacuum  conditions  demon¬ 
strated  heat  transfer  coefficients  closer  to  the  Nusselt 
theory  than  when  only  five  active  tubes  per  tube  bundle 
column  were  considered. 
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The  results  of  all  runs  yielded  the  following  exponent 
for  the  following  Nusselt  type  equations: 


r*% 

—  =  Nr  -  (N 

h. 
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0  .  82 

61 

tt 

0.23 

0.75 

62 
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0.21 

0.78 

for  a  tube-to-tube  spacing  of  3  mm,  condensate  drops 
barely  detached  from  the  upper  tube  before  striking 
the  next  lower  tube  at  low  condensation  rates .  At 
higher  condensation  rates,  drops  elongated  to  the 
next  lower  tube,  instantaneously  linking  the  tubes. 

For  both  the  above  cases ,  the  condensate  drops  were 
rapidly  integrated  into  the  condensate  film  of  the 
lower  tube  rather  than  rolling  over  the  tube  perimeter 
With  increasing  rates  of  condensation,  the  number  of 
drops  per  unit  length  of  a  given  tube  increased.  When 
the  number  of  drops  per  unit  length  on  the  top  tube 
was  compared  to  the  number  of  drops  per  unit  length 


on  each  of  the  lower  tubes,  it  was  found  that  the 
ratios  remained  constant  despite  changes  in  the  rate 
of  condensation. 
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VII.  RECOMMENDATIONS 

A.  TEST  APPARATUS  MODIFICATIONS 

To  improve  the  uncertainty  and  validity  of  the  test  results, 
the  following  test  apparatus  modifications  are  considered 
advisable : 

1.  The  most  significant  problem  which  remains  with  the 
test  condenser  is  the  breakage  of  tube  wall  thermo¬ 
couples.  Due  to  the  cost  of  the  thermocouples  and 
the  problems  involved  in  replacing  damaged  ones,  it 
would  be  best  to  use  larger  diameter,  larger  wall 
thickness  tubes .  This  would  allow  for  the  in-house 
fabrication  of  thermocouples  from  cheaper  fiberglass 
insulated  thermocouple  wire,  and  for  the  installation 
of  thermocouples  to  the  tube  grooves  using  epoxy  cement 
and  copper  strips  rather  than  having  to  crimp  and 
silver  solder  then  in  place  as  was  the  current  practice. 

2.  The  test  condenser  should  be  properly  leveled  and  piping 
targeted  to  the  test  condenser  to  reduce  or  eliminate 
the  condenser  tube  leveling  problem. 

3.  The  test  condenser  flanges  should  be  made  wider  to  allow 
for  more  sealing  area  which  would  eliminate  the  air 
leakage  problem  along  the  flanges. 

4.  For  the  Nusselt  curves,  the  N=5  position,  for  both  nor¬ 
malized  heat  transfer  coefficients,  occurs  near  the 

6  6 


knee  of  the  curve.  It  is  considered  that  more  reliable 
data  would  be  obtained  if  more  tubes  were  added  to 
the  tube  bundle.  Three  additional  tube  rows  of  16  mm 
O.T.  tubes  could  be  installed  in  the  existing  rig. 

It  would  be  feasible,  with  some  piping  modi f ication , 
to  further  increase  the  number  of  tube  rows  by  adding 
a  flanged  section  directly  below  the  existing  conden¬ 
ser  section.  The  existing  steam  supply  may  limit  the 
number  of  tube  rows  that  can  be  installed  and  still 
allow  for  a  sufficiently  large  volume  and  velocity  of 
exhaust  steam.  It  is  also  considered  advisable  to 
install  a  row  of  dummy  tubes  beneath  the  last  active 
tube  for  use  when  additional  tubes  are  to  be  simulated. 

5 .  The  existing  cooling  water  rotameters  should  be  replaced 
with  lower  maximum  flow  rate  rotameters.  This  would 
allow  for  the  attainment  of  lower  flow  rate?  .  a  lower 
cooling  water  AT,  and  improved  uncertainty  of  the 
results.  All  runs  in  this  study  were  made  at  about 

15%  maximum  flow  rate  (corresponding  to  a  cooling 
water  velocity  of  about  1.5  m/s)  which  was  about  the 
lowest  setting  one  could  use  with  good  accuracy. 

6.  For  simulation  of  larger  tube  bundles,  the  inundation 
water  should  be  at  the  temperature  of  the  condensate 
falling  from  the  bottom  tube  of  the  test  condenser. 
Therefore  ,  instrumentation  should  be  installed  to 
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to  measure  this  temperature.  Additionally,  the  tem¬ 
perature  control  system  for  the  porous  tube  water 
supply  tank  should  be  modified  to  allow  for  more  rapid 
heating  and  cooling  to  avoid  long  delays  between  run 
increments . 

?.  Installation  of  a  window  on  the  back  of  the  test  con¬ 
denser  would  allow  for  viewing  of  the  back  side  of 
the  tubes  so  that  the  absence  of  drop-wise  condensa¬ 
tion  could  be  continually  verified. 

ADDITIONAL  TESTS 

The  following  additional  tests  should  be  conducted, 

1.  To  assure  proper  removal  of  non-condensable  gases  during 
atmospheric  runs,  the  condenser  exhaust  steam  was  dumped 
to  the  bilges  thus  preventing  vapor  velocity  calculations. 
A  set  of  runs  should  be  made  with  the  exhaust  steam 
piped  to  the  secondary  condenser,  the  vent  on  the  second¬ 
ary  condenser  open,  and  the  air  ejector  secured.  The 
requirement  for  the  securing  of  the  air  ejector  is 
necessitated  since  it  can  pull  a  vacuum  on  the  test 
condenser  by  itself,  however,  if  a  larger  vent  wet? 
installed  on  the  secondary  condenser  hotwell,  it  may 
be  possible  to  operate  with  the  air  ejector  in  opera¬ 
tion.  The  results  of  the  above  tests  should  then  be 
compared  to  the  data  presented  herein.  If  build-up 
of  non-condensable  gases  can  be  avoided,  then  all  steam 


could  be  condensed  and  the  vapor  velocity  calculated. 
Alternatively,  a  venturi  type  velocity  measuring 
instrument  could  be  installed  in  the  steam  piping 
between  the  desuperheater  and  the  test  condenser. 

2 .  Runs  should  be  made  varying  the  cooling  water  flow 
rates  so  that  a  comparison  could  be  made  between  the 
method  used  in  this  study  and  the  Wilson  Plot  technique. 

3.  Verification  of  the  results  of  this  study  should  be  con¬ 
ducted  by  making  sample  runs  with  a  full  set  of  tube 
wall  thermocouples  installed  in  each  active  tube. 

4.  Movies  should  be  made  of  the  condensation  process 
process  so  that  accurate  conclusions  could  be  drawn 

as  to  condensate  drop  phenomena  and  their  relationship 
to  condenser  performance. 

5 .  Conduct  additional  tests  to  ascertain  the  effect  of 
the  change  in  vapor  velocity  through  the  tube  bundle 
on  the  heat  transfer  coefficients. 

6 .  Conduct  tests  with  enhanced  tubes  after  a  good  smooth 
tube  performance  base-line  has  been  verified. 
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Table  I . 


Channel  Numbers  for  Stainless  Steel  Sheathed 
Copper-Constantan  Thermocouples 
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Comparison  of  Heat  Transfer  Coefficients 
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Table  VIII. 


Comparison  of  Heat  Transfer  Coefficients 
and  Mean  Vapor  Velocity  for  a  Single 
Uninundated  Tube  Under  Vacuum  Conditions 
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Figure  3:  Mean  Flow  Width  for  a  Square  In-Line  Tube  Bundle 
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AREA  OF  UNIT  CELL  =  /I  (?Td)  / 2 

FREE  AREA  OF  CELL  =  d2  ( /I  Prp2/2  -  7T/4) 

MEAN  FLOW  WIDTH  =  d  ( ?T  -  * / 2  /3  PT) 

(  =  FREE  AREA/HEIGHT  OF  CELL) 


ure  4:  Mean  Flow  Width  for  a  Staggered  Tube  3undle 


TUBE  NUMBER 

COMPARISON  OF  PREDICTIONS  OF  TUBE  BUNDLE  PERFORMANCE 


Figure  6  :  COMPARISON  OF  PREDICTIONS  FOR  THE  LOCAL  HEAT  TRANSFER  COEFFICIENT 
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igure  Jb:  Rear View  of  Test  Apparatus 


KEY  TO  FIGURES  7a  AND  7b 


1.  Test  Condenser 

2 .  Test  Condenser  Hotwell 

3 .  Secondary  Condenser 

4.  Secondary  Condenser  Hotwell 

5 .  Test  Condenser  Steam  Exhaust  Piping 

(Small  pipe  connected  for  atmospheric  pressure  run. 
Larger  pipe  is  used  for  vacuum  runs.) 

6.  Porous  Tube  Supply  Piping  (top  pipe)  and  Cooling  Water 
Supply  Piping  (lower  five;  note  the  cooling  water  inlet 
thermocouple  connections.) 

7.  Cooling  Water  Rotameters 

3.  Porous  Tube  Supply  Water  Tank 

9 .  Porous  Tube  Supply  Water  Pump 

10 .  Porous  Tube  Rotameter 

11.  Porous  Tube  Supply  Water  Temperature  Controller 

12.  Test  Condenser  Manometer 

13.  Desuperheater  Tank  with  Water  Spray  Nozzles 

14 .  Desuperheater  Water  Supply  Tank 

15 .  Desuperheater  Water  Supply  Pump 
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NOTE !  ALL  COMPONENTS  DRAWN  TO  SCALE 
Figure  3:  Sketch,  of  Test  Condenser 
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gure  9:  Detailsof  Transition  Piece  and  Vortex  Annihilator 


.  EXHAUST  (STEAM  AND  CONDENSATE) 
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STEAM  TO  SECONDARY 
CONDENSER  HOTWELL 


1 


TO  TEST  CONDENSER  HOTWELL 


gure  10:  Details  of  Exhaust  and  Condensate  Pining  from  the 
Exhaust  Plenum 
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16  mm  DUMMY  T 


16  mm  ACTIVE  T 


Figure  13:  Schematic  Diagram  of  Cooling  Water  and  Porous  Tube  Water  Systems. 
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Figure  14 :  Schematic  Diagram  of  Condensate  and  Vocuum  Systems. 
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Figure  15:  Schematic  Diagram  of  Steam  System. 


Figure  16:  Location  as  Thermocouples. 
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Figure  20:  AVERAGE  HN/H1  VS  TUBE 


Figure  21:  AVERAGE  HN/H1  VS  TUBE 


Figure  22:  AVERAGE  HN/H1  VS  TUBE  NUMBER 


Figure  23:  AVERAGE  HN/Hl  VS  TUBE 


Figure  24:  AVERAGE  HN/Hl  VS  TUBE  NUMBER  RUNS  21 


AVERAGE  HN/H1  VS  TUBE  NUMBER  RUNS  31 
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Figure  28:  AVERAGE  HN/H1  VS  TUBE  NUMBER  RUNS  41 


Figure  30  :  AVERAGE  htl/Hl  VS 


Figure  31:  AVERAGE  HN/H1  VS  TUBE  NUMBER 


Figure  33:  AVERAGE  W/Hl  VS  TUBE  NUMB*  BUS  51 


Figur  e  34:  AVERAGE  HN/H1  VS  TUBE 
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Figure  36  :  HN/H1  VS  TUBE 


Figure  37:  HN/HI  VS  TUBE  MJMBER  RUN  61 


Figure  38;  HN/H1  VS  TUBE 


Figure  39:  HN/H1  VS  TUBE 


APPENDIX  A:  OPERATING  PROCEDURES 

INITIAL  PROCEDURES 

1.  Energise  the  main  circuit  breaker  located  in  Power 
Panel  P-2  on  the  wall  to  the  right  of  the  test 
apparatus . 

2.  Energize  the  circuit  breaker  on  the  left  side  of  the 
old  control  board  by  pressing  the  ON  button. 

3.  Energize  the  following  switches  in  the  control  panel: 

a.  HI  -  Porous  Tube  Water  Supply  Pump  (Feed  Pump) 

b.  #2  -  Outlets 

c.  H 3  -  Hot  Water  Heater 

d.  #4  -  Condensate  Pump 

e.  # 6  -  Cooling  Tower 

f.  HI  -  Cool  Water  Pump 

4.  Ensure  all  test  apparatus  valves  are  closed. 

5 .  Fill  the  desuperheater  and  porous  tube  condensate 
supply  tanks  with  distilled  water.  Set  the  tempera¬ 
ture  controller  for  the  porous  tube  water  supply  tank 
at  about  50°C,  fully  open  the  recirculation  valve, 
P-1,  and  start  the  pump  to  begin  heating  the  water. 
The  controller  will  have  to  be  reset  to  the  proper 
supply  temperature  once  steady  state  conditions  are 
obtained . 
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6.  Fill  the  cooling  water  supply  tank.  This  can  be  done 
by  the  use  of  a  hose  connected  to  the  house  water 
supply,  by  backfilling  with  valves  CK1-1  and  CW-4 
open,  or  by  opening  valve  WS-1  and  the  water  supply 
valve  on  the  steam  hot  water  heater. 

7.  Start  the  air  compressor  following  the  posted 
instructions . 

3.  Energize  the  instrumentation  as  per  Appendix  3. 

OPERATION 

1.  House  Steam 

a.  Open  the  main  supply  valve. 

b.  Open  valve  MS-3  until  the  pressure  gage  indicates 
the  desired  steam  supply  pressure. 

c.  Fully  open  MSD-1  to  drain  any  condensate. 

d.  Open  valves  MS-4  and  MS-5  until  the  desired 
steam  supply  pressure  is  obtained  and  readjust 
MS-3  as  necessary. 

2.  Condensate  and  Vacuum  Systems 

a.  Open  valve  C-4  to  supply  cooling  water  to  the 
secondary  condenser. 

b.  Open  valve  A-l  to  energize  air  ejector. 

c.  To  collect  the  condensate  in  the  test  condenser 
hotwell,  operate  the  system  with  valve  C-l 
closed.  After  a  test  run  is  completed,  open 
valve  C-l  to  drain  the  condensate  into  the 
secondary  condenser  hotwell. 
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d.  Whenever  the  secondary  condenser  hotwell  level 
approaches  half  full,  the  hotwell  must  be  pumped 
a^»jn  to  maintain  vacuum.  To  pump  down  the  hotwell, 
open  valves  C-2  and  C-3  and  start  the  condensate 
pump.  Secure  the  pump  down  by  securing  the  con¬ 
densate  pump  and  closing  valves  C-2  and  C-3. 

3.  Cooling  Water  System 

a.  Open  valves  CW  -  j. ,  CW-2  ,  and  CW-3. 

b.  Ensure  valves  CK1-1  and  CW-4  are  closed. 

c.  Energize  the  two  cooling  water  pumps. 

d.  Open  valves  CW-5 ,  CW-6 ,  CW-7,  CW-8 ,  and  CW-9  to 
obtain  the  desired  cooling  water  rates. 

4.  Desuperheater  System 

a.  Start  the  supply  tank  pump. 

b.  Open  the  recirculation  valve  D-l  one  turn. 

c.  Set  the  desuperheater  rotameter  to  a  10-20%  flow 
rate  by  adjusting  valve  D-3. 

d.  Fully  open  valve  D-4 . 

e.  Depending  upon  the  house  steam  state  point,  the 
desuperheater  may  or  may  not  be  required.  The 
optimum  test  condenser  steam  inlet  temperature 
is  such  that  the  degree  of  superheat  is  less 
than  10°C.  The  system  is  not  required  when  run¬ 
ning  the  test  apparatus  under  atmospheric  pres¬ 
sure  conditions.  The  temperature  is  adjusted  by 
regulating  the  flow  rate  of  water  to  the  desuper¬ 
heater  spray  nozzles  via  the  rotameter. 
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5 .  Porous  Tube  Water  Supply  System 

a.  Once  steady  state  conditions  have  been  achieved 
for  a  run,  reset  the  temperature  controller  to 
the  proper  inundation  temperature  . 

b.  Adjust  the  rotameter  to  the  required  flow  rate 
for  each  run.  The  supply  tank  recirculation  valve 
may  have  to  be  adjusted  to  achieve  the  desired 
flow  rate,  but  should  never  be  fully  closed  since 
pump  damage  and  non-uniform  water  temperature 

may  occur. 

c.  Refill  the  supply  tank  as  required  with  distilled 
water  to  maintain  the  tank  level  above  the  heating 
element . 

6 .  Miscellaneous 

To  maintain  a  clear  test  condenser  window,  open  valve 
A-2  and  then  energize  and  adjust  the  air  heater  power 
supply.  When  securing,  always  turn  off  the  power 
supply  first  and  allow  the  air  heater  to  cool  before 
securing  valve  A-Z. 

C.  SECURING  THE  TEST  APPARATUS 

1.  Secure  the  steam  valves  MS-5,  MS-4,  MS-3,  and  the 
main  supply  valve  . 

2.  Secure  the  air  compressor. 

3.  Secure  the  desuperheater  by  securing  the  pump  and 
valves  D-l  and  D-3.  Drain  the  supply  tank  by  opening 
valve  D-2  (if  desired). 
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Secure  the  porous  tube  water  supply  system  by  securing 
the  pump,  temperature  controller,  and  valves  P-1  and 
P-4.  Drain  the  supply  tank  by  opening  valve  P-3  (if 
desired) . 

Secure  the  air  ejector  by  closing  valve 

Secure  the  test  condenser  viewing  window  air  heater 

as  prescribed  above. 

Secure  the  instrumentation. 

Allow  the  test  section  to  cool  down  for  about  one-half 
hour,  then  secure  the  coding  water  pumps  and  close 
valves  C3-1,  CW-2  ,  CW-3,  CW-5  ,  CW-6  ,  CW-7,  CW-8,  CW-9, 
and  C-4 . 

Drain  the  test  condenser  and  secondary  condenser 
hotwells  . 

Secure  all  circuit  breakers . 

Drain  the  cooling  water  system  piping  and  rotameters 
by  opening  valve  CW-10  and  cracking  open  the  cooling 
water  rotameter  supply  valves. 

Drain  the  cooling  water  supply  tank  by  opening  the 
drain  valve  via  the  remote  operating  rod. 


APPENDIX  3:  AUTODATA  MINE  SCANNER  OPERATION 


1.  Plug  the  power  cord  into  the  wall  outlet. 

2.  Turn  the  main  power  switch  on.  The  pressure  transducer 
power  supply  has  a  separate  switch,  located  inside  the 
back  door  of  the  Autodata  Nine  Scanner,  which  must  also 
be  turned  on. 

3.  Set  the  time. 

a.  Place  all  alarms  and  output  switches  in  the  off 
position . 

b.  Set  the  data  and  time  on  the  thumbwheels. 

c.  Set  the  display  switch  to  "time.” 

d.  Lift  the  "set  time"  switch. 

4.  Assigning  Multiple  Channels 

a.  Set  the  display  switch  to  "all." 

b.  Check  that  all  alarms  and  switches  are  "off." 

c.  Set  the  scan  switch  to  "continuous." 

d.  Lift  the  "slow"  switch. 

e.  Set  the  first  channel  thumbwheels  to  "000"  and  the 
last  channel  thumbwheels  to  "001." 

f.  To  assign  channel  "000"  and  "001,"  depress  and  hold 
the  "10V"  and  "HI  RES"  buttons  for  at  least  one  scan 
and  lift  the  scan  start  switch  to  start  scanning. 
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g.  Set  the  last  channel  thumbwheels  to  ’'0  39  "  setting 
the  first  channel  thumbwheels  to  ’’002  ." 

h.  Depress  the  "skip”  button  and  lift  the  scan  switch 
to  skip  channels  00Z  through  039. 

i.  Set  the  last  channel  thumbwheels  to  "039"  before 
setting  the  first  channel  thumbwheels  to  "040." 

j  .  To  assign  channels,  depress  and  hold  the  "T/°C"  and 
"HI  RES”  buttons  for  at  least  one  computer  scan. 

k.  Set  the  scan  switch  to  the  fast  scan  position. 

5 .  Interval  Scan 

a.  Set  the  thumbwheels  to  the  time  interval  desired 
between  scans  . 

b.  Depress  the  "stop  enter"  switch. 

c.  Set  the  display  switch  to  "interval." 

d.  Depress  the  "set  interval"  switch. 

e.  Set  the  scan  switch  to  "interval." 

f.  Set  the  first  channel  thumbwheels  to  "000." 

g.  Set  the  last  channel  thumbwheels  to  "089." 

h.  Lift  the  "scan  start"  switch. 

6  .  Optional  Devices  as  Required 

a.  Printer  on/off. 

b.  Single  channel  display. 
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APPENDIX  C:  UNCERTAINTY  ANALYSIS 


The  general  form  of  the  Kline  and  McClintock  [Ref.  23] 
"second  order"  equation  is  used  to  calculate  the  uncertainty. 
If  the  resultant,  R,  is  some  function  of  the  primary  variables 
X,  X^ ,  •••,  Xn,  then  the  uncertainty  is  given  by. 


(C-l) 


where  SX^,  5X2 ,  ...,  6Xn  are  the  uncertainties  in  each  of  the 

measured  variables  X-^,  X2,  .  .  .  ,  Xn . 

1.  Uncertainty  in  the  Steam  Side  Heat  Transfer  Coefficient, 
h  from  Equation  (11): 


h  = 


m  C  (T  -  T  .  ) 
p  co  ci 

A  (Ts  -  Tw> 


Applying  Equation  (C-l)  to  the  above,  the  following 
equation  is  obtained: 
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The  following  uncertainties  were  assigned  to  the  variables: 

6m  =  0.01  .kg/s 

<5C  =  0.004  kj/kg  •  °C 

6T  =  0.1  °C 

c° 

6T  . =  o.r  °c 

Cl 

5A  =  0.0001  m2 

5Ts  =  0.5  °C 

6T~  =  1.0  °C 

W 

2.  Uncertainty  in  the  Normalized  Average  Heat  Transfer 
Coefficient,  Kj’j'/h^ 

The  normalized  average  heat  transfer  coefficient  is  obtained 
for  the  Nth  tube  by  taking  the  average  of  the  heat  transfer 
coefficients  of  the  first  N  tubes  and  dividing  this  by  the  heat 
transfer  coefficient  of  the  first  tube: 

h^  (h^  +  hj  +  ...  +  h^)/N 


Applying  Equation  (C-l)  to  the  above,  the  following 
equation  results: 


1/2 


where 

NT  =  the  tube  number. 

The  uncertainties  for  5h^  and  Sh^  (i  =  1,  2,  . .  .N)  are  obtained 
from  Equation  (C-2).  For  example: 


5<K?V  ; 

[(  5hi  \2 

(  Shz  ) 

‘V'V 

L\h1+h2+h3+h4+h5  / 

'hl+1"12+h3+h4+h5 

(  4h3  )2  * 

(  ) 

\hi+h2+h3+h4+h5 

'  h-^+h2  +h 3+h4+hg  / 

(  6hs  -  Y  * 

1/2 

(6Jh)2l 

^  Vh2+h3+h4+hS  ; 

V  J 
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3.  Uncertainty  in  the  Averaged  Normalized  Average  Heat 

Transfer  Coefficient,  .. 

(hN/h1) 

This  coefficient  is  simply  h^/h^  averaged  for  X  number 
of  runs : 


(  -  )  =  i  i  (  hJL  ) 

hn  X  i=l  \  ,  ' 


‘1  i 


Applying  Equation  (C-l)  to  the  above,  the  following 
equation  results: 


6(hN/h1) 


1/2 

rx  I 

Li=l  'z  <Vhl>i  ’  ■* 


i  =  l 


for  N  >  2 


for  N  =  1 


(C-4) 


The  uncertainties  for  the  <5(h^/h^)^  terms  are  obtained  from 
Equation  (C-3).  For  example  for  the  third  tube  and  five  runs 


N  =  3  and  X  =  5 


« (TTJ/h1) 

(hj/h^ ) 


7  *<Wl 

A  nwi 


/fiCKT/h.  ),■  \2 

(rwi) 


+ 


/  &(K3/hih  V 

\  FwT  / 


/  6(h3/hl)4  V 

v  rw  / 
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1/2 


«VVs_\! 

*'  £  (KI/hl>i 

'i  =  l 


where  E  (h,/h, ) .  = 
i  =  l  1  1 


CA-  & 


*1 


1 '3 


C  C  —  5  ) 


©. 


*7 


1  '5 


4.  Uncertainties  for  the  Normalized  Local  Heat 
Transfer  Coefficient,  h^/h^ 

This  coefficient  is  simply  the  heat  transfer  coefficient 
of  a  given  tube,  N,  divided  by  that  of  the  first  tube  or, 
for  the  fifth  tube,  N=S  and: 

hN  .  h5 
HI  "  hi 

An  application  of  Equation  (C-l)  results  in  the  following 
equation : 

1/2 

for  N  >  2 


5(hN/h1) 


(W 


r 

2  *  (Sh»U  21 

< 

-  hi 

{^>  J 

( C-6 ) 


for  N  =  1 


Where  the  5h^  and  Sh^  values  are  obtained  from  Equation  (C-2). 
For  example: 
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6(h2/h1) 

\rl i\2  ♦ 

(Sh2')2l 

(h^/h^) 

Ll  hx  ’ 

Uj  J 

1/2 


5(h5/h1) 

r/Shi  \2  *  /5hs  \2 1 

<W 

L\  hi  /  U5  /  ■* 

5 .  Uncertainties  for  the  Averaged  Normalized  Local 
Heat  Transfer  Coefficient,  (n^77TTp 

This  coefficient  is  simply  the  h^/h^  values  averaged  for 

X  number  of  runs: 


Applying  Equation  (C-l)  to  the  above,  results  in  the 
following  equation  (which  is  similar  to  Equation  (C-4): 

<5  (hM/h,  ) 

N  1 

(H^/h^) 

0  for  N  =  1 


1/2 


=  < 


for  N  >  2 


( C-  7 ) 


Where  the  6(h^/h^)^  values  are  obtained  from  Equation  (C-6). 
Note  that  Equation  (C-7)  has  the  exact  form  as  Equation  (C-4) 
but  with  the  (h^/h-^K  terms  replaced  by  (h^/h^)^  terms. 
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5.  Uncertainties  in  Musselt's  Single  Tube  Heat 
Transfer  Coefficient,  h^u 

From  Equation  (1): 


h..  =  .72  5 

Nu 


[ 


-  -  °v>  hfi;  8 

U  D  <TS  -  Tw) 


r 


Applying  Equation  (C-l)  to  the  above  and  considering  Pv 
negligible  in  relation  to  P,  the  following  equation  results: 


The  following  uncertainties  were  assigned  to  the  variables 

ok  =  ±0.001  w / m  •  °C 

5p  =±0.1  kg/m2 

5hftj  =  ±0.1  kJ/kg 

5g  =  ±0.001  m/s2 

5u  =  ±  1  x  10'6  kg/m-s 

5D  =  ±0.0001  m 
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5T  =  ±1.0°C 
s 

5T  =  ±1.0°C 
w 


7 .  Uncertaint 
From  Equation 


h 

o 


10 . 74 


ies  in  Fujii's  Correlation 
(9)  one  can  obtain: 


where 


U 


00 


7 

m 


m 


COND  v 


mf 


from  Equation  (12). 


Applying  Equation  (C-l)  to  the  equation  for  V  one  obtains 
the  following  equation: 


/{"C0MD  \  2  ,  /f^_¥  + 

\  mC  0  MD  J  \  v  / 


( C-9 ) 


The  following  uncertainties  were  assigned  to  the  variables: 

5mC0MD  =  ±0-0005  kS/s 
5v  =  ±0.01  m3/kg 

5A  -  =  ±0 .0001  m2 

mf 

Applying  Equation  (C-l)  to  Equation  (9),  as  modified 
above,  results  in  the  equation: 
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Applying  Equation  ( C  —  1 )  to  the  above  results  in  the  equation: 


or 


1/2 

<Sq  /  A  _  ^  /  <Sq  V  +  /  5 A 

q/A  L\  q  /  \  A 

where 

<5q  is  obtained  from  Equation  (C-ll) 
6A  =  0  .0001  m2 . 


( C-12 ) 


Tables  C-I  and  C-II  present  uncertainty  values  for  various 
runs.  From  the  tables,  the  following  uncertainties  can  be 
assigned  to  the  calculated  results  of  this  study: 

For  atmospheric  pressure  runs: 

5h  =  ±0.7  kw/m2  •  °C 

=  ±0.07 
hl 

6r-±i  =  ±0.07 
hl 

t0-02 


16  U 


±0.02 


r 


•II 

*  - 

°o 

±0.1  kw/m2  • 

°C 

5q/A  i 

±20  kw / m2 

vacuum 

runs  : 

oh 

=  ±2.3  kw/m*- 

•  °C 

hl 

=  ±0.17 

•S 

i  ±0.17 

i  ±0.07 

•C?) 

=  ±0.03 

5h.. 

Mu 

=  ±0.5  kw/m2 

.  oc 

6h 

o 

i  ±0.3  kw/m2 

•  °C 

<5V 

m 

=  ±0.4  m/s 

<5q/A 

=  ±10  kw/m2 
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Uncertainties  for  Runs 
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t'1,  i  iM~mriai«r*i 11 ^  i 


Table  C-II 


Selected 

Uncertainties 

for  Top 

Tubes 

±5hNu 

±6h 

o 

±  5V_ 
m 

±6  q  /  A 

RUN  # 

(kw/m2  *°C) 

(kw/m2  *°C) 

(m/s) 

(kw/m2 ) 

1 

0  .1 

-- 

— 

20 

11 

0.1 

-- 

-- 

20 

21 

0.1 

— 

— 

20 

31 

0.1 

-- 

20 

41 

0.5 

0.5 

0  .  4 

9 

46 

0  .5 

0 . 3 

0  .  3 

9 

51 

0.5 

0.3 

0.4 

10 

5  6 

0.5 

0.3 

0.3 

10 
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APPENDIX  D:  SAMPLE  CALCULATIONS 


The  following  illustrates  the  calculation  procedures  of 
this  study.  Tubes  1  and  2  of  runs  1  and  46  were  chosen  as 
examples.  The  data  used  came  from  Tables  II  and  III;  while, 
the  uncertainties  were  obtained  from  Appendix  C  and  Tables 
C-I  and  C-II  thereto. 

A.  ATMOSPHERIC  PRESSURE  RUN  EXAMPLE,  RUN  1 

1.  q  =  m  C  ( T-T  • ) ,  q/A  =  m  C  ( T  -T  • ) /A 

n  p  CO  Cl  ^  p  CO  Cl 

-  Tube  1 

m  =  10.977  kg/min  =  0.18295  kg/s 
T  =  3  3  . 5°C  T  •  =  2  7 . 1°C 

CO  Cl 

TS  =  (Tco+Tci)/2  =  30'3°C 
Cp  (at  30  . 3°C)  =  4.130  kJ/kg-°C 

q  =  4.8942794  kw 

A  =  tt (  .  3048)  (  .0159  )  m2 

q/A  =  321.45983  kw/m2 

<5q/A  =  t20  kw/m2 

q/A  =  320  ±  20  kw/m2 

Tube  2 

m  =  0.18295  kg/s 

Tco  =  33*°°C  Tci  =  27.2°C  Tb=30.1°C 

Cp  =  4.180  kJ/kg-°C 
q  =  4.4354398  kw 
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CK  =  7T(  .  30  4  8 )  (  .0159  )  m 
q/A  =  291.32297  kw/m2 
6q/A  =  ±20  kw/ra2 
q/A  =  290  ±  20  kw/m2 

h  =  (Ts-fw) 

Tube  1 

T  =  101°C  T  =  7  8 . 2°C 

s  w 

2 

q/A  =  320  ±  20  kw/m 
h1  =  14.035083  kw/m2 -°C 
6hq  =  1.1  kw/m2  *°C 

=  14.0  ±  1.1  kw/m2 -°C 
Tube  2 

T  =  101°C  T  =  7 0  . 7°C 

s  w 

q/A  =  290  ±  20  kw/m2 
h2  =  9.5709571  kw/m2  *°C 
<5h2  =  0.7  kw/m2  -°C 
h2  =  9.6  ±  0.7  kw/m2 -°C 

^N^l  =  ^22^1 

h2/h1  =  9.6/14.0  =  0.6857 

6h2/h1  =  ±0.07 

h2/h1  =  0 .68  ±  0 .07 

S/hl  Z  ^2*^1 

KJ  =  (hq  +  h2)/2  =  (14.U  9 . 6  )  / 2  =  11.8  kw/m 

hj/h1  =  11.8/14.0  =  0.8428 
6h^’/h1  =  ±0.08 

h^/h1  =  0.84  ±  0.08 
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(h^/h-^)  for  runs  1- 

10 ,  tube  2 

RUN  h2/h1 

RUN 

Vhi 

1  0.68 

6 

0  .77 

2  0.62 

7 

0.67 

3  0.71 

8 

0  .69 

4  0.76 

9 

0.68 

5  0.73 

10 

0  .69 

(Data  for  runs  2-10 

taken 

from 

Table  V) 

10 

(h2/h1)  =  j-g  i  =  1 

(h2/h 

l>i  = 

7.05 

10 

6h2/h]_)  =  ±0.02 

(h2/h1)  =  0.70  ± 

0.02 

(H^/h^)  =  (h'J/h^) 

RUN  hj/l^ 

RUN 

1  0.84 

6 

0  .38 

2  0.31 

7 

0  .  84 

3  0.86 

8 

0.85 

4  0.88 

9 

0  .  84 

5  0189 

10 

0  .  85 

(Data  for  runs  2-10 

taken 

from 

Table  V) 

10 

(h2/h1)  =  ro  i=l 

(*?hl 

)  .  = 
i 

8.54  . 

Hnr  '  0 

5(ITJ/h1)  =  ±0.02 


(H^/h1)  =  0.85  ±  0.02 
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0.705 


.  854 


f 


1/4 


h.,  for  tube  1 
Nu 


hMu  =  -725 


9  .725 


k  P(0-Pv)  hfg 


“  D  (Ts  -  Tw 


1/4 


V 


3  2  , 

P  hfg  g 


y  D  (Ts  -  Tw) 


T  =  10 1°C  T  =  7  8 . 2°C 

s  w 

=  (T  +  T  )/2  =  89 . 6°C 
r  s  w 

2 

k  =  0.674  w/m • °C  p  =  967.6  kg/m 

hf  =  2,254,300  J/kg  p  =  3.19  x  10-4  kg/m-s 

g  =  9.91  m/s2  D  =  0 .0159  m 

h„  =  11.094  kw/m2 -°C 

6h„  =  ±0.1  kw/m2  *°C 

Nu 

h„  =  11.1  ±  0 .1  kw/m2 -°C 
Nu 


B.  VACUUM  RUN  EXAMPLE ,  RUN  46 

1.  q  =  *  Cp(Tco  -  Tci),  q/A  =  m  Cp(T00  -  Tci>/A 
Tube  1 

m  =  0.18295  kg/s 

T  o  =  2  2  . 2°C  T  •  =  19 . 7°C  T0  =  21.0°C 

Cp  =  4.182  kJ /kg ■ °C 

q  =  1.9127423  kw 

A  =  tt (  .  3048)  (  .0159)  m2 

q/A  =  125.53033  kw/m2 

5q/A  =  ±10  kw/m2 

q/A  =  130  ±  10  kw/m2 
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Tube  2 


m  =  0.13295  kg/s 
T  =  2  2 . 4°C  T  .  =  19 . 7°C 

CO  Cl 

Cp  =  4.132  kJ/kg  *°C 
q  =  2  .0657616  kw 
A  =  tt(  .  30  4  8 )  (  .0159  )  m2 
q/A  =  135 .68075  kw/m2 
5q/A  =  ±10  kw/m2 

q/A  =  140  ±  10  kw/m2 

2.  h  =  q/A  (Ts  -  Tw) 

Tube  1 

T  =  5  2  . 2°C  T,  =  44 . 9°C 

s  w 

q/A  =  130  ±  10  kw/m2 
h1  =  17 . 30  8219  kw/m2  *°C 
6h1  =  ±3.0  kw/m2 -°C 
hx  =  18.0  ±  3.0  kw/m2  -°C 
Tube  2 

TS  =  52*2°C  ^  =  42 . 2°C 

q/A  =  140  ±  10  kw/m2 
h2  =  14.0  kw/m2 -°C 
5h2  =±2.1  kw/m2 -°C 
h2  =  14.0  ±  2.1  kw/m2 -°C 

3.  hN/hl  =  h2/h1 

h2/h1  =  14.0/18.0  =  0.77777778 
6h2/h1  =  ±0.17 

h2/h1  =  0.78  ±  0.17 
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. 


(13.0  +  14 . 0 ) / 2 


4  .  h^/h-^  =  /h^ 

hj  =  (h1  +  h2)/2  = 

h^/hL  =  16.0/18.0  =  0.88888889 

5h2/h1  =  ±0.13 

h^/h1  =  0.89  ±  0.13 

c .  (h^/h^)  for  runs  46-50,  tube  2 


RUN 

h2 

/h 

46 

0  . 

78 

47 

0  . 

80 

48 

0  . 

32 

49 

0  . 

84 

50 

0  . 

88 

=  16  . 


(Data  for  runs  47-50  taken  from  Table  V) 


_  5 

(h2/h1)  =  j  if1  (h2/h1)i  =  4.12/5'=  0.824 


6(h2/h1)  =  ±0.08 


(h2/h1)  =  0.82  ±  0.08 


6.  (FT^/h^)  for  runs  45-50  ,  tube  2 


RUN 

h. 

/h 

46 

0 

39 

47 

0 

90 

43 

0 

.91 

49 

0 

.92 

50 

0 

.34 

(Data  for  runs  47-50  taken  from  Table  V) 

_  5 

(H^/hL)  =  j  il1  (E^/h1)i  =  4.56/5  =  0. 

6(h^/h1)  =  ±0.09 

( In  /h,  )  =  0  .91  ±  0  .09 

W  X 


0  kw/rn^ • ° 


912 
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7. 


for  tube  1 


3  2 

K  P 


hfg-  -5-- 1 

(T  -  T  )J 


=  52.2  C 


44 ,9°C 


=  48.6  C 


k  =  3.633  w/m-  C 


p  =  990.7  kg/ m" 


hf  =  2,377,500  J/kg  u  =  5.76  x  10~  kg/m-s 

g  =  9.81  m/s2  D  =  0.0159  m 

hv.  =  12.519  kw/m2  -°C 
Nu 

5h,r  =  ±0.5  kw/m2*°C 
Nu 

h„  =  12 .5  ±  0 .5  kw/m2 -°C 
Nu 


Vm  =  mC0ND  V/Amf 

Amf  =  Nc  d  (PT  -  */4V  L 


N  =31/3 
c 

=  1.5 


d  =  0.0159  m 


L  =  0.3048  m 


A  f  =  0.01577318  m 
6Amf  =  ±0.0001  m2 

A  f  =  0 .0158  ±  0  .0001  m2 

raC0ND  =  +  ^SC;  x  10  ^  P 


=  4.32  ml/s  =  7.69  ml/s 

p  =  990.7  kg/m2 

=  [(4.32  s-  7.69)  x  10-6  m3/s]  990.7  kg/m3 


‘COND 


=  1.1876  x  10"  kg/s 


6lflC0ND  =  ±0  •0005  k§/s 


0.0119  ±  0.0005  kg/s 


/ 


V 


m 
<5  V 


m 


V. 


m 


9  .  h 


10.947  m  /kg 
3.24489  m/s 
±0.3  m/s 
8.2  ±  0.3  m/s 


k. 


=  10.74 


/V  d  \ 


0 . 312 


k. 

i-i 

Vm 

ho 

6h, 


0 .638  w/m-  C 


8.2  ±  0.3  m/s 
20.135  kw/m2  *°C 
±0  .  3  kw/m2 -°C 


=  0.0159  m 
o 

vL  =  0.581  x  10“15  m2/s 


20.1  ±  0.3  kw/m 


2  Or 
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APPENDIX  E:  COMPUTER  PROGRAM  AND  DOCUMENTATION 

The  computer  program  listed  in  Table  E-l  uses  the  HP-85 
computing  system  and  basic  language.  The  program  uses  an 
inter-active  approach. 

A.  DIRECTIONS/QUESTIONS  OF  THE  COMPUTER  PROGRAM 

The  following  is  a  list  of  directions/questions  and 
responses  required  of  the  operation: 

1.  Line  50  -  "Enter  run  number."  The  opertor  inputs  the 
number  designation  of  the  run. 

2 .  Line  80  -  "What  is  the  number  of  tubes?"  The  operator 
inputs  the  number  of  active  tubes  simulated  in  the 
experiment . 

3.  Line  130  -  "Enter  water  and  steam  temperatures." 

The  operator  inputs  for  each  tube  the  two  cooling 
water  inlet  temperatures,  the  three  cooling  water 
outlet  temperatures,  and  the  steam  saturation 
temperature . 

4.  Line  150  -  "Enter  tube  wall  temperatures."  The 
operator  inputs  for  each  tube  the  top,  bottom,  and 
averaged  center  wall  temperatures . 

5.  Line  170  -  "Enter  cooling  water  mass  flow  rate  and 
specific  heat."  The  operator  inputs  the  cooling 
water  mass  flow  rate  and  specific  heat  for  each 
tube  . 
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Items  3-5  will  be  displayed  for  each  active  tube. 

3.  COMPUTATIONS  OF  THE  PROGRAM 

1.  Lines  210-240  -  The  average  cooling  water  inlet  and 
outlet  temperatures  are  calculated. 

2.  Lines  310-315  -  The  heat  transferred  and  heat  flux 

are  calculated  by: 

q  =  m  C  (T__  -  T  .)  and, 

H  p  co  ci  ’ 

q/A  =  m  C  (T  -  T  . )/A 
n  p  CO  Cl 

3.  Line  330  -  The  outside  neat  transfer  coefficient  is 
calculated  by: 

h  =  - a - 

A  'Ts  -  V 

4.  Lines  340-350  -  The  average  heat  transfer  coefficient 
is  calculated. 

5.  Line  360  -  The  normalized  heat  transfer  coefficient 
is  calculated. 


C. 


DATA 

.  OUTPUT 

1. 

Line 

70 

2  . 

Line 

250 

3  . 

Line 

250 

4  . 

Line 

270 

5  . 

Line 

320 

6  . 

Line 

370 

7  . 

Line 

380 

3. 

Line 

381 

Run  number. 

-  Average  cooling  water  outlet  temperature . 

-  Average  cooling  water  inlet  temperature. 

-  Average  tube  wall  temperature . 

-  Tube  number. 

-  Normalized  heat  transfer  coefficient. 

-  Heat  transferred. 

-  Heat  flux. 


130 


Line  382 


The  outside  heat  transfer  coefficient. 


9  . 

10.  Line  383  -  The  average  heat  transfer  coefficient. 

D.  GRAPHICS  OUTPUT 

The  graphics  are  plotted  using  a  peripheral  plotter. 

The  operator  must  ensure  that  a  3  1/2"  x  11"  sheet  of  paper 
is  in  the  plotter.  The  program  will  stop  before  plotting 
the  graphics  and  will  not  commence  the  graphics  output  until 
the  "continue"  key  is  pressed. 
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TABLE  E-l :  COMPUTER  PROGRAM 

10  1  CONDI  PLOTS  HN/H1  VRS  TUBE  NUMBER 

20  OPTION  3ASE  1 

30  DIM  Q( 30 ) ,  M( 30 ) ,  11(30),  12(30),  01(30),  02(30), 

03(  30),  1(  30),  0(  30),  K( 30 )  ,  Tl(30),  T2(30),  T3(30) 
40  DIM  H ( 3  0 ) ,  T0 ( 30 )  ,  T4(30),  T5(30),  T6(30),  T( 30)  , 

C(  30)  ,  Q 1 ( 30) 

50  DISP  "ENTER  RUN  NUMBER" 

60  INPUT  V 

70  PRINT  "RUN  NUMBER  IS,"  V 

80  DISP  "WHAT  IS  THE  NUMBER  OF  TUBES?" 

90  INPUT  J 

100  A1  =  PI  *  .3048  *  .0159 

110  K1  =  0 

120  FOR  N  =  1  TO  J 

130  DISP  "ENTER  WATER  AND  STEAM  TEMPERATURES" 

140  INPUT  Il( N)  ,  I2(N),  01 (N) ,  02(N),  03(N),  T0(N) 

150  DISP  "ENTER  TUBE  WALL  TEMPERATURES" 

160  INPUT  T4(N),  T5(N),  T6(N) 

170  DISP  "ENTER  COOLING  WATER  MASS  FLOW  RATE  AND 

SPECIFIC  HEAT" 

180  INPUT  M(N)  ,  C(N) 

190  NEXT  N 
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200 

FOR  N 

=  1  TO  J 

2  10 

0  ( N )  = 

( 01 ( N )  +  0  2 ( N )  +  03(N) 

220 

I(M)  = 

(I1(N)  +  I2(N))/2 

230 

T(N)  = 

(T4(N)  +  T5(N)  +  T6(N) 

250 

PRINT 

"0  =  ",  0(N) 

2  6  0 

PRINT 

"I  =  ",  I(M) 

270 

PRINT 

"T  =  ",  TOO 

280 

PRINT 

"T0  =  ”,  T0(N) 

290 

NEXT  N 

300 

FOR  N 

=  1  TO  J 

310 

Q(N)  = 

M(N)  *  C(N)  *  ( 0 (N )  - 

315  ' 

Q1(N) 

=  Q(N)/A1 

320 

PRINT 

"TUBE  NUMBER  =  ",  N 

330 

K(N)  = 

Q(N)/A1/(T0(N)  -  T( N) ) 

340 

K1  =  K1  +  K(N) 

350 

K2  =  Kl/N 

360 

H(N)  = 

K2/KC1) 

370 

PRINT 

"N.H.T.C.  =  ",  H(N) 

380 

PRINT 

"Q  =  ",  QCN) 

381 

PRINT 

"Q/A  =  ",  Q1(N) 

332 

PRINT 

"H  =  ",  K(N) 

383 

PRINT 

"HH(AVG)  =  ",  K2 

385 

NEXT  N 

390 

DISP  " 

ENSURE  PAPER  IN  PLOTTER 

400 

PAUSE 

183 


410 


PLOTTER  IS  705 


420  LIMIT  32,  235,  28,  182 

430  LOCATE  10,  130,  IS,  96 

440  SCALE  1,  30,  0,  1.1 

450  FXD  0,1 

460  LAXES  1,  .1,  1,  0,  2,  .1 

470  SETGU 

430  MOVE  71,  7.5 

490  DEG  @  LDIR  0 

500  LORG  5 

510  LABEL  "TUBE  NUMBER" 

520  MOVE  3,  56.5 

530  DEG  @  LDIR  90 

540  LORG  5 

550  LABEL  "HN/H1" 

560  MOVE  1,  55 

570  LDIR  90 

530  LORG  5 

590  LABEL  " _ " 

600  MOVE  71,  2.5 

610  LDIR  0 

620  LORG  5 

630  LABEL  "FIG.  HN/H1  VS  TUBE  NUMBER  RUN  #,"  V 

640  MOVE  58,  4.5 

6  50 


LDIR 


660 


LORG  5 


670 

LABEL  " _ " 

580 

SETUV 

690 

MOVE  1,  1 

700 

FOR  T  =  1  TO  30  STEP 

.1 

710 

R  =  T  ~  -  .25 

720 

DRAW  T,  R 

730 

NEXT  T 

740 

MOVE  1,  1 

750 

FOR  T  =  1  TO  30  STEP 

.1 

760 

R  =  .6+  . 42  *  T  ~  - 

.25 

770 

DRAW  T,  R 

730 

NEXT  T 

790 

FOR  N  =  1  TO  J 

800 

MOVE  N,  HCN) 

810 

LABEL 

820 

NEXT  N 

330 

END 

185 


LIST  OF  REFERENCES 


1.  Standards  for  Steam  Surface  Condensers,  6th  ed.,  Heat 
Exchange  Institute,  1 3 7 0  . 

2 .  Standards  of  Tubular  Exchanger  Manifacturers  Association, 
4th  ed . ,  Tubular  exchanger  Manif acturers  Association, 

Inc .  ,  1959  . 

3.  Search,  H.  T.,  A  Feasibility  Study  of  Heat  Transfer 
Improvement  in  Marine  Steam  Condensers,  MSME ,  Naval 
Postgraduate  School ,  Monterey ,  California ,  December  1977. 

4.  Beck,  A.  C.,  A  Test  Facility  to  Measure  Heat  Transfer 
Performance  of  Advanced  Condenser  Tubes,  MSME ,  Naval 
Postgraduate"  School,  Monterey,  California ,  January  19  77  . 

5.  Pence,  D.  T.,  An  Experimental  Study  of  Steam  Condensation 
on  a  Single  Horizontal  Tube,  MSm£ ,  Naval  Postgraduate 
School,  Monterey,  California,  March  1978. 

6.  Reilly,  D.  J.,  An  Experimental  Investigation  of  Enhanced 
Heat  Transfer  on  Horizontal  Condenser  Tubes,  MSME,  Naval 
Postgraduate  School,  Monterey,  California,  March  1978. 

7.  Fenner,  J.  H.,  An  Experimental  Comparison  of  Enhanced 
Heat  Transfer  Condenser  Tubing,  MSME,  Naval  Postgraduate 
School , ~ Monterey ,  California ,  September  1978. 

3.  Ciftci,  H.,  An  Experimental  Study  of  Filmwise  Condensation 
on  Horizontal  Enhanced  Condenser  Tubing,  MSME ,  Naval 
Postgraduate  School ,  Monterey ,  California ,  December  1979. 

9.  Marto ,  P.  J.,  Reilly,  D.  J.,  Fenner,  J.  H.,  "An  Experi¬ 
mental  Comparison  of  Enhanced  Heat  Transfer  Tubing," 
Advances  in  Enhanced  Heat  Transfer,  Proceedings  of  the 
13th  National  Conference,  American  Society  of  Mechanical 
Engineers,  New  York,  New  York,  1979. 

10.  Demire  1,  I.,  The  Effects  of  Condensate  Inundation  on 
Condensation  Heat  Transfer  in  lube  Bundles  of  Marine 
Condensers,  MSMli,  Naval  Postgraduate'  School,  Monterey , 
California ,  December  1980. 

11.  Eshleman,  D.  E. ,  An  Experimental  Investigation  of  the 
Effect  of  Condensate  Inundation  on  Heat  Transfer  in.  "a 
Horizontal  "lube  Bundle,  M5Me!  ,  Naval  Postgraudate  School , 
Monterey,  California,  March  1980. 


136 


•ftU- 


* 


12.  Morrison,  R.  H.,  A  Test  Condenser  to  Measure  Condensate 
Inundation  Effects~in  a  Tube  Bundle,  MSME,  Naval  Post¬ 
graduate  School,  Monterey,  California ,  March  1981. 

13.  Nobbs,  D.  W.,  The  Effect  of  Downward  Vapour  Velocity  and 
Inundation  on  the  Condensate  Fates  on  Horizontal  Tubes 
and  Tube  Fank~  Ph.D  Thesis,  University  of  Bristol, 
Bristol,  England,  April  19  75  . 

14.  Eissenberg,  D.  M.,  An  Investigation  of  the  Variables 
Affecting  Steam  Condensation  on  the  Outside  of  a  Horizon¬ 
tal  Tube  Bundle ,  Ph. D  Thesis  ,  University  of  Tennessee, 
Knoxville  ,  Tennessee ,  December  1972. 

15.  Jakob,  Max,  Heat  Transfer,  Vol.  1,  J.  Wiley  S  Sons,  Inc., 
1948. 

16.  Kutateladze,  S.  S.,  Heat  Transfer  in  Condensation  and 
Boiling,  Translated  from  a  publication  of  the  State 
Scientific  and  Technical  Publishers  of  Literature  on 
Machinery,  Moscow-Leningrad ,  1952  for  the  United  States 
Atomic  Energy  Commission. 

17.  Shklover,  G.  G.  and  Buevich,  A.  V.,  "The  Mechanism  of 
Film  Flow  with  Steam  Condensation  in  Horizontal  Bundles 
of  Tubes,"  Teploenergetika ,  25  (  4)  62-65  ,  1978  . 

18.  Yung,  D.,  Lorenz,  J.  J.,  and  Ganic,  E.  N.,  "Vapor/Liquid 
Interaction  and  Entrainment  in  Falling  Film  Evaporators," 
Transactions  of  ASME,  Vol.  102,  February  19 80. 

19.  Fuks,  S.  H.,  "Heat  Transfer  with  Condensation  of  Steam 
Flowing  in  a  Horizontal  Tube  Bundle,"  Teploenergetika , 
4(1)  35-39,  1957.  English  Translation  NEL  1041, 

National  Engineering  Laboratory,  East  Kilbride. 

20.  Chen,  M.  M.,  "An  Analytical  Study  of  Laminar  Film  Conden¬ 
sation:  Part  2  -  Single  and  Multiple  Horizontal  Tubes," 

Journal  of  Heat  Transfer,  February  1965. 

21.  Fujii,  T.,  "Vapor  Shear  and  Condensate  Inundation," 
Research  Institute  of  Industrial  Science,  Kijushu 
University,  Fukuoka,  Japan,  1979. 

22.  Knudsen,  J.  G.  and  Katz,  D.  L.  ,  Fluid  Dynamics  and  Heat 
Transfer,  pp.  332-335,  McGraw-Hill  Book  Company,  1958 . 

23.  Kline,  S.  J.  and  McClintock,  F.  A.,  Describing  Uncer¬ 
tainties  in  Single  Sample  Experiments,  Mech.  Engin., 

Vol .  74,  pp .  3-8,  January  1953. 


187 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Department  Chairman,  Code  69  2 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  p.  J.  Marto,  Code  69Mx  5 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Professor  R.  H.  Nunn,  Code  69Nn  1 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 

6.  Mr.  R.  W.  Kornbau  2 

Code  2721 

David  Taylor  Naval  Ship  Research  and 
Development  Center 
Bethesda,  Maryland  20034 

7.  LT  P.  J.  Noftz ,  USN  (N411)  2 

Commander  Naval  Surface  Forces, 

U.S.  Atlantic  Fleet 
Norfolk,  Virginia  23511 

3.  Commander  1 

Naval  Sea  Systems  Command 
Navy  Department 
Washington,  DC  20362 


183 


